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PREFACE 
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Validation  Study  of  the  CRREL  Shallow  Snow  Mobility  Model 


GEORGE  L.  BLAISDELU  PAUL  W.  RICHMOND,  SALLY  A.  SHOOP, 
CHARLES  E.  GREEN  AND  RUSSELL G.  ALGER 


INTRODUCTION 

In  1986,  CRREL  developed  a  model  for  predict¬ 
ing  the  mobility  of  tracked  and  wheeled  vehicles 
operating  on  shallow  snow  covers.  (Shallow  snow 
is  defined  here  by  requiring  that  the  deformational 
zone  in  the  snow  beneath  the  tire  or  track  intersect 
the  ground  or  firm  surface  upon  which  the  snow 
lies.  In  contrast,  a  deep  snow  appears  to  a  vehicle  to 
be  semi-infinite.)  This  shallow  snow  mobility  model 
predicts  the  basic  mobility  parameters  of  motion 
resistance  and  gross  driving  traction.  From  these 
quantities,  net  traction,  or  the  amount  of  tractive 
force  available  for  acceleration,  payload,  towing 
force,  hill  climbing,  etc.,  can  be  calculated. 

Tire  initial  version  of  the  shallow  snow  model 
used  an  empirically  derived  relationship  between 
snow  shear  strength  and  initial  density  to  make  its 
predictions  of  gross  traction.  Tlris  relationship  has 
been  developed  from  mobility  data  obtained  with 
the  CRREL  Instrumented  Vehicle  (CIV)  over  past 
field  seasons.  Motion  resistance  is  predicted  from  a 
calculation  of  the  energy  expended  in  compacting 
the  snow  using  a  theoretically  established  expres¬ 
sion.  During  January  1988,  a  2-year  mobility  study 
was  launched  jointly  by  WES  and  CRREL.  Tlris  so- 
called  zviitter  zvheels/tracks  study  provided  the  first 
opportunity  to  compare  the  model's  predictions 
with  actual  snow  mobility  data  for  a  wide  variety  of 
vehicles.  It  was  also  the  intent  of  this  study  to  com¬ 
pare  directly  the  performance  of  wheeled  and 
tracked  vehicles  in  snow .  Tlie  selected  family  of  test 
vehicles  included  four  wheeled  vehicles  and  three 
tracked  vehicles. 


During  the  zolieels/ tracks  study,  we  also  hoped  to 
add  to  the  CIV  snow  strength^ensity  data  base, 
and  to  explore  whether  or  not  the  strength-density 
relationship  developed  with  the  CIV  could  be  du¬ 
plicated  using  a  shear  annulus  device.  Tliis  report 
describes  the  results  of  our  first  year's  effort. 

BACKGROUND 

The  primary  goal  of  this  study  was  to  validate 
the  shallow  snow  mobility  model.  As  previously 
mentioned,  this  model  is  based  on  theoretical  rela¬ 
tionships  and  empirical  expressions  developed  in 
the  past,  and  a  large,  but  scattered,  data  base.  Since 
the  shallow  snow  model  has  not  appeared  in  the 
literature  in  the  past,  we  begin  by  describing  the 
approach  used  in  version  1 .0  (SSMl  .0)  for  predict¬ 
ing  mobility. 

The  two  principal  quantities  that  govern  mobil¬ 
ity  are  traction  and  motion  resistance.  Commonly, 
two  types  of  traction  (net  and  gross)  and  two  types 
of  motion  resistance  (internal  and  external)  are 
identified. 

External  motion  resistance  R^,  for  our  case,  is  a 
manifestation  of  the  energy  expended  in  deform¬ 
ing  and  displacing  the  snow.  Internal  motion  resis¬ 
tance  is  attributable  to  a  combination  of  tire  flexing 
or  friction  on  track  linkages,  friction  on  the  bear¬ 
ings  and  gears  along  the  vehicle's  drive  line,  and 
other  inefficiencies  inherent  in  transferring  torque 
from  the  engine  to  the  snow  surface  (or  other  op¬ 
erating  surface).  Most  resistance  measuring  tech¬ 
niques  unavoidably  result  in  a  measure  of  the  total 


motion  resistance  (both  external  and  internal). 
Our  interest  is  in  understanding  and  predicting 
only  the  external  resistance  due  to  snow  defor¬ 
mation,  since  internal  resistance  has  no  fundamen¬ 
tal  significance  in  limiting  traction. 

Gross  traction  T  is  defined  as  the  forward  thrust 
developed  at  the  ^re/ snow  or  track/snow  inter¬ 
face.  It  is  a  function  of  the  snow's  shear  strength 
and  the  geometric  and  adhesive  ability  of  the  tire  or 
track  to  engage  this  shear  strength.  The  tractive 
device's  ability  to  pull  a  load,  accelerate,  or  climb  a 
slope  is  found  by  summing  the  forces  at  the  contact 
patch.  Thus 


fn  ~  “  Rs 


(1) 


where  is  the  net  traction  and  is  a  measure  of  the 
traction  available  for  use  by  the  vehicle.  Here,  T  is 
positive  pointing  towards  the  direction  of  intenc^d 
travel  and  R^is  positive  in  the  opposite  direction. 

All  traction  measurement  methods  known  to  us 
measure  T^.  These  tests  are  also  designed  to  deter¬ 
mine  the  maximum  net  traction  for  a  given 
tire-snow  or  track-snow  condition,  and  it  is  this 
maximum  value  that  is  most  often  reported.  If  we 
assume  that  the  vast  majority  of  tires  or  tracks  in 
use  are  capable  of  engaging  the  shear  strength  of 
the  snow  (usually  true  for  all  but  bald  tires  or 
smooth  tracks),  then  it  follows  that  the  snow's 
shear  strength  places  an  upper  bound  on  T  .  From 
eq  1  then,  the  same  is  true  for  T^.  Thus,  a  reasonable 
approach  for  predicting  maximum  is  to  base  it 
on  the  shear  strength  of  snow. 

Using  energetics  principles,  we  can  describe  the 
compaction  of  snow  by  a  tire  or  track.  Equating  the 
energy  of  compaction  to  resistance  to  vehicle  mo¬ 
tion  in  snow,  we  can  obtain  a  theoretical  expres¬ 
sion  for  R  . 

S 

The  shallow  snow  model  SSMl.O  described  in 
Appendix  A  uses  this  approach  to  produce  esti¬ 
mates  of  T  and  R  and  to  calculate  net  traction  (i.e., 
mobility)  ^om  eq  1 . 


FIELD  EXPERIMENTS 

Test  location  and  test  sites 

The  field  experiments  described  in  this  report 
were  conducted  at  the  Keweenaw  Research  Center 
(KRC)  located  at  the  Houghton  County  Airpark, 
Michigan.  KRC  is  located  on  the  Keweenaw  Penin¬ 
sula  of  upper  Michigan,  about  900  km  (550  miles) 
northwest  of  Detroit.  Two  vehicle  test  areas  were 
used  to  obtain  differing  snow  conditions.  The 


majority  of  the  testing  took  place  at  the  KRC  "ice 
rink,"  a  large  (25,000-m^[270,000-ft^])  leveled  area 
often  used  for  winter  vehicle  traction  tests  on  ice. 
The  remaining  tests  were  done  on  an  unplowed 
taxiway  at  the  Houghton  County  airport. 

Test  vehicles 

CRREL  imtrumeuted  vehicle 

The  CRREL  instrumented  vehicle  (CIV)  is  based 
on  a  1977  Jeep  Cherokee  with  modifications  to  its 
braking  and  driving  components  to  accommodate 
typical  mobility  tests.  In  addition  to  an  on-board 
computer-based  data  acquisition  system,  the  ve¬ 
hicle  contains  a  number  of  transducers  for  force 
and  speed  measurements.  The  front  two  wheels  are 
equipped  with  triaxial  load  cells  that  measure  the 
vertical,  longitudinal  and  lateral  forces  acting  at  the 
tire/ground  interface.  Individual  wheel  and  ve¬ 
hicle  speeds  can  also  be  measured.  For  this  test 
program,  the  CIV  was  equipped  with  Goodyear 
Tiempo  tires,  size  P225/75R1 5,  and  they  were  tested 
with  inflation  pressures  of  179  and  103  kPa.  Basic 
vehicle  and  tire  characteristics  are  given  in  Appen¬ 
dix  B. 

Wheels/tracks  test  vehicles 

The  vehicles  and  test  tires  used  in  the  winter 
wheels/ tracks  performance  program  were: 

M988  High-Mobility  Multipurpose  Wheeled 
Vehicle  (HMMWV),  4x4,  MicheUn  37x1 2.5R1 6.5LT 
tires. 

LAV  25  Light  Armored  Vehicle,  8x8,  Michelin 
11.00R16  tires. 

LAV  25  Light  Armored  Vehicle,  8x8,  Michelin 
12.50R20  tires. 

M977  Heavy  Expanded  Mobility  Tactical  Truck 
(HEMTT),  8x8,  Michelin  16.0R20  tires. 

M923  5-ton  truck,  6x6,  Goodyear  14.00R20XL 
tires.* 

M973  Small  Unit  Support  Vehicle  (SUSV),  ar¬ 
ticulated,  tracked. 

M113A1  Armored  Personnel  Carrier  (APC), 
tracked. 

M2  Bradley  Fighting  Vehicle  System  (BFVS), 
tracked. 

The  gross  vehicle  test  weights  and  individual  tire 
or  track  loads  are  given  in  Appendix  B. 


*  The  M923  was  not  initially  selected  for  the  ivheels/tracks 
program;  however,  it  was  used  because  of  its  availabil¬ 
ity,  and  because  it  differed  from  the  other  wheeled 
vehicles.  It  was  not  loaded  to  its  published  capacity. 
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Test  procedures 

The  test  procedures  followed  are  typical  for 
mobility  field  studies,  in  that  measurements  of  net 
traction  T^and  total  motion  resistance  (internal 
plus  external  resistance)  were  made  with  each 
vehicle  under  varying  snow  conditions.  Wheeled 
vehicle  tests  were  generally  conducted  at  two  tire 
inflation  pressures  to  document  the  effect  of  tire 
deflection  on  mobility  levels.  A  shear  annulus  device 
was  also  operated  side  by  side  with  the  vehicles. 
Snow  type,  depth,  temperature  and  density  were 
measured  as  well  ciS  air  temperature  and  solar 
input. 

CIV 

Motion  resistance  is  measured  with  the  CIV 
while  the  rear  tires  are  driving  and  the  front  wheels 
are  rolling  free.  Since  the  triaxial  load  cells  are 
located  just  inside  the  front  wheels,  this  test  meas¬ 
ures  the  amount  of  resistance  felt  at  the  front  tires. 
With  the  rear  tires  driving,  the  whole  vehicle  mass, 
including  the  front  axles  from  the  center  of  the 
vehicle  out  to  the  load  cells,  is  being  moved  by  the 
rear  wheels.  Force  measurements  are,  therefore, 
against  this  frame  of  reference  and  constitute  total 
front  tire  motion  resistance  R^only. 

The  internal  motion  resistance  R.  is  first  estab¬ 
lished  on  a  level,  undeformable  surface.  Measure¬ 
ments  of  are  obtained  for  each  tire  type  and 
selected  inflation  pressure.  By  convention,  motion 
resistance  tests  are  done  at  a  vehicle  speed  of  8  km/ 
hr.  Our  past  test  experience  in  snow  has  shown  that 
resistance  values  are  independent  of  slight  vari¬ 
ations  in  vehicle  speed.  Variations  in  values 
between  tire  types  are  the  result  of  differences  in 
the  forces  necessary  for  tire  flexing  and  can  be 
attributed  to  differences  in  their  design. 

Motion  resistance  tests  on  snow,  then,  result  in 
measurement  of  both  internal  and  external  motion 
resistance  R^.  R^is  found  by  subtracting  from  the 
value  of  Rj  measured  when  the  CIV  is  operating  on 
snow. 

A  traction  test  is  performed  with  the  CIV  bv 
accelerating  tl  front  (driving)  wheels  while  brak¬ 
ing  the  rear  wheels  to  hold  the  vehicle  speed  con¬ 
stant  at  8  ±  1  km/hr  (5  ±  0.5  miles/hr).  In  this 
manner  the  front  tires  are  driven  through  a  wide 
range  of  slip  values,  starting  at  zero  slip.  A  plot  of 
measured  T^vs  Differential  Interface  Velocity  (DIV) 
is  used  to  obtain  tractive  effort.  The  value  is 
taken  as  an  average  of  the  upper  1 5%  of  the  tractive 
force  readings.  This  results  in  averaging  over  a 
range  of  slip  values  that  could  reasonably  be  main¬ 
tained  by  a  vehicle  operator.  Gross  traction  is 


then  calculated  using  eq  1  for  each  pair  of  resistance 
and  traction  tests. 

Wheels j tracks  vehicles 

Tractive  performance  of  the  wheels/tracks  ve¬ 
hicles  was  determined  by  measuring  DrawBar  Pull 
(DBF)  and  Towed  Motion  Resistance  (TMR).  Draw¬ 
bar  pull  tests  were  conducted  by  measuring  the 
force  that  a  vehicle  can  exert  on  a  cable  that  is  being 
used  to  resist  vehicle  motion.  Thus 

DBF  =  (Tn)  W)  (2) 

where  N  is  the  number  of  driven  tires  or  tracks. 
Motion  resistance  is  determined  by  measuring  the 
vehicle's  resistance  to  towing  with  an  instrumented 
cable.  Here 

TMR  =  (R,)  (h)  (3) 

where  n  is  the  number  of  tires  or  tracks. 

Our  procedure  for  measuring  DBF  is  as  follows. 
A  load  vehicle  of  approximately  the  same  size  and 
performance  as  the  test  vehicle  was  selected  to 
apply  a  resistance  to  the  test  vehicle.  A  steel  cable 
0.016  m  (5/8  in.)  in  diameter  and  15  m  (50  ft)  in 
length  was  connected  from  the  front  of  the  load 
vehicle  to  a  load  cell  attached  to  the  rear  pintle  hook 
of  the  test  vehicle.  A  string  pay-out  system  for 
measurement  of  true  ground  distance  was  also 
mounted  on  the  test  vehicle.  Tachometers  mounted 
on  the  drive  wheels  or  sprockets  were  used  to 
measure  wheel  or  track  travel  during  a  test. 

Each  test  vehicle  was  operated  in  its  lowest  gear 
and  at  optimum  engine  rpm  (yielding  a  ground 
speed  of  3-8  km/hr  [2-5  miles/hr]).  The  vehicle 
proceeded  into  the  test  lane  with  the  load  vehicle 
following  and  the  cable  between  the  two  vehicles 
being  slack  and  unloaded.  The  load  vehicle  driver 
initially  applied  load  to  the  test  vehicle  by  braking 
and  then  increased  the  applied  load  in  steps.  The 
test  sequence  proceeded  from  the  test  vehicle  ini¬ 
tially  experiencing  a  no  load-no  slip  condition  and 
increased  by  steps  up  to  a  high  load-high  slip 
condition  or  a  power  limited  condition  in  which  the 
test  vehicle  could  not  maintain  the  desired  track  or 
tire  speed.  Forward  speed  and  wheel  or  track  speed 
were  maintained  (very  nearly)  constant  for  several 
seconds  of  steady-state  pull  measurements.  Ve¬ 
hicle  slip  was  calculated  from  a  magnetic  tape 
record  of  bo  th  vehicle  and  ground  speed .  Tlie  vehicle 
slip,  calculated  for  short  time  windows  in  the  rec¬ 
ord,  in  percent,  is  equal  to 
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Slip  =  (100).  (4) 

rfv 

where  =  apparent  distance  traveled  by  the 
wheel  or  track 

=  distance  actually  traveled  by  the  ve¬ 
hicle. 

Continuous  measurements  were  made  in  this 
manner  until  a  sufficient  number  of  load  and  slip 
combinations  were  recorded  to  develop  a  drawbar 
pull-slip  curve  (usually  two  complete  test  se¬ 
quences). 

The  procedure  used  for  obtaining  TMR  of  the 
test  vehicle  in  each  test  area  was  to  tow  the  vehicle 
(with  its  transmission  in  neutral)  at  a  speed  of 
approximately  3  km/hr  (2  miles/hr).  After  each 
traction  test,  the  vehicle  was  steered  into  an  undis¬ 
turbed  area  adjacent  to  each  traction  test  area, 
usually  in  a  position  straddling  the  ruts  of  the 
associated  traction  test.  The  load  vehicle  then  towed 
the  test  vehicle  to  determine  TMR.  The  test  pro¬ 
ceeded  for  a  sufficient  distance  to  permit  the  load 
cell  readings  to  stabilize  and  to  be  recorded  on 
magnetic  tape.  An  average  value  during  the  stable 
portion  of  the  record  was  taken  as  TMR.  TMR  was 
also  measured  on  an  undeformable  surface.  R  is 

S 

once  again  found  by  taking  the  difference  between 
measurements  on  snow  and  the  hard  surface. 

As  with  in  the  CFV  traction  test,  DBF  is  a 
function  of  slip  or  DTV.  Generally,  in  low  density 
snow  (less  than  500  kg/m^)  maximum  DBP occurs 


at  rates  of  slip  greater  than  20%.  However,  little 
useful  work  is  being  done  at  very  high  slip  rates 
and  the  efficiency  of  operation  is  inversely  propor¬ 
tional  to  slip.  Thus,  the  maximum  DBP,  which  oc¬ 
curred  near  20%  shp,  was  used  for  the  calculation 
of  T^(eq  2).  Equation  1,  again,  is  used  to  determine 
gross  traction,  where  R^  is  obtained  from  the  tow¬ 
ing  tests. 

Shear  annulus  device 

The  KRC  shear  annulus  used  in  this  study  is 
installed  on  a  slide  rail  that  is  mounted  on  the  side 
of  a  Pole  Cat.  This  allows  a  series  of  tests  to  be  run 
without  moving  the  vehicle.  A  computer  program 
(run  on  an  IBM  PC)  controls  the  start  and  stop  of 
data  collection  from  outside  the  vehicle,  and  also 
collects  data  throughout  the  sequence  of  tests. 

The  annulus  used  in  this  study  had  a  smooth 
rubber  face.  The  outside  radius  of  the  ring  was  9.21 
cm  (3.625  in.)  and  the  inside  radius  was  6.67  cm 
(2.625  in.).  The  surface  area  of  the  rubber  ring  is 
then  126.64  cm^  (19.63  in.^).  Figure  1  shows  the 
shear  device. 

The  test  procedure  involves  rotating  the  annu¬ 
lus  at  10  rpm,  while  recording  rotational  displace¬ 
ment  and  applied  torque.  Each  set  of  shear  tests  is 
run  using  five  different  normal  loads  on  the  annu¬ 
lus.  These  loads  are  in  increments  of  89  N  (20  Ibf) 
and  range  from  89-445  N  (20-100  Ibf).  The  first  run 
is  made  without  any  added  weight,  since  the  shaft 
and  annulus  weigh  89  N.  The  apparatus  is  lowered 
to  the  surface  of  the  snow  and  the  initial  sinkage  of 


Figure  1.  Field  shear  annulus  device. 
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the  annulus  is  measured.  An  outside  operator  then 
controls  the  start  and  stop  of  the  test  from  a  remote 
switch.  Once  the  test  has  started,  the  annulus  is 
allowed  to  turn  through  one  full  revolution.  After 
one  revolution,  the  operator  releases  the  switch 
and  both  the  motor  and  the  computer  data-acc]ui- 
sition  system  stop.  The  final  sinkage  of  the  ring  is 
measured  and  the  annulus  is  pulled  out  of  the 
snow.  After  each  run,  the  device  is  moved  forward 
on  the  slide  rail  to  a  new  position.  Another  89-N 
increment  of  weight  is  added  and  the  test  is  re¬ 
peated.  When  all  five  runs  are  complete,  the  data 
arestored.  The  entire  sequence  of  five  normal  loads 
is  repeated  three  times  at  each  test  site. 

RESULTS 

CIV  traction  and  motion  resistance 

Idealized  traction  test  results  for  a  driven  tire  on 
the  CIV  are  ^hown  in  Figure  2.  The  CIV  directly 
measures  during  a  traction  test.  The  average 
value  of  the  top  15%  of  the  data  (to  remove  local 
spikes  in  the  force  record)  is  taken  as  the  peak 


Shear  t./isplacement  or 
Rate  of  Shear  Displacement 


Figure  2.  Idealized  traction  test  results. 

tractive  force.  A  typical  resistance  test  result  is 
shown  in  Figure  3,  where  R^is  taken  as  an  average 
of  all  of  the  data. 

Table  1  summarizes  the  data  obtained  with  the 
CIV  during  the  ivheels/tracks  field  program  at  KRC. 
Date,  temperature,  snow  density  and  depth,  and 


Table  1.  CIV  traction  and  motion  resistance  data. 


Dnte 

Tei’ifi. 

Density 

(kg/w3) 

Depth 

(cm) 

Pressure 

(kPa) 

Measured 

gross 

traction 

coefficient 

TglW 

Measured 

external 

resistance 

coefficient 

RJW 

19  Jan  88 

-2 

550 

5 

179 

0,46 

0.020 

103 

0.43 

0.032 

21  Jan  88 

-5 

60 

2-5 

170 

0.36 

0.014 

103 

0.35 

0.023 

22  Jan  88 

-5 

160 

7 

170 

0,38 

0.005 

103 

0,37 

0.019 

23  Jan  88 

-5 

75 

12 

170 

0.33 

0.019 

103 

0.35 

0.028 

23  Jan  88 

-2 

75 

10 

170 

0.32 

0.010 

103 

0.33 

0.023 

26  Jan  88 

-5 

’  :s 

9 

170 

0.35 

0.048 

103 

0,34 

0.056 

10  Feb  88 

-8 

220 

31 

170 

0.37 

0.164 

103 

0.41 

0.192 

10  Feb  88 

-7 

250 

22 

170 

_ 

0.116 

103 

0.37 

0.129 

12  Feb  88 

-16 

230 

22 

170 

0.36 

0.115 

103 

0.45 

0.161 

5 


Distance  Traveled 


Figures.  Idealized  resistance  test  results. 

tire  pressure  are  self  explanatory  (a  consensus  of 
values  was  taken  when  KRC  and  CRREL  records 
differed).  R  and  T  (for  two  driven  wheels)  are 
listed  in  coefficient ^orm  (R^/ Wand  T  /W).  Each 
traction  and  resistance  coefficient  reported  in  Table 
1  is  the  average  of  at  least  two  traction  tests  (in  most 
cases  four  to  six  tests).  The  CIV's  average  normal 
wheel  load  is  6174  N  (1388  Ibf). 

Wheels! tracks  vehicles  traction 
and  motion  resistance 

Table  2  summarizes  the  traction  and  motion 
resistance  data  taken  for  the  wheeled  and  tracked 
vehicles  used  in  this  study.  The  values  listed  for 
motion  resistance  and  peak  gross  traction  are  in 
coefficient  form  (R^/ W  and  7"^/ W,  respectively), 
and  were  obtained  from  the  rneasured  TMR  and 


DB  P  values  (eq  2  and  3,  see  Appendix  C  for  sample 
test  vehicle  field  results).  Resistance  and  traction 
from  the  test  vehicle  records,  although  for  whole 
vehicles  rather  than  for  single  wheels,  are  deter¬ 
mined  in  the  same  manner  as  was  done  for  the  CIV 
(see  Fig.  2  and  3),  and  will  be  discussed  later. 

Shear  annulus  device 

Figure  4  shows  how  measured  shear  stress  var¬ 
ies  with  time  during  a  shear  annulus  device  test. 
From  these  data,  the  shear  strength  is  assigned 
the  value  of  the  constant  stress  level  x  that  results 
after  an  initial  spike  caused  by  static  friction.  The 
plot  of  shear  strength  values  for  the  various  normal 
loads  shown  in  Figure  4  is  given  in  Figure  5.  The 
slope  of  a  least  squares  regression  line  through 
these  data  is  defined  as  tan(l)  and  the  y-intercept  of 
this  plot  is  called  cohesion  c.  For  the  test  shown,  a 
straightlinefits  thedata  very  well.  This  is  represen¬ 
tative  of  most  of  the  tests  run.  Table  3  is  a  complete 
listing  of  the  results  over  the  test  period  (see  Table 
4  for  associated  snow  properties). 

Accuracy  and  limitations  of  data 

The  CIV's  load  cells  are  accurate  to  ±  2%  of 
measured  value  down  to  130  N  (30  Ibf).  Forces 
lower  than  this  are  shown  from  calibration  to  be 
accurate  to  ±  9  N  (2  Ibf). 

DBP  and  TMR  of  the  wheels/ tracks  vehicles  were 
measured  using  load  cells  with  the  following  ranges: 
BFVS,  HEMTT— 89,000  N  (20,000  Ibf) 

LAV,  SUSV,  AFC,  5-ton-44,500  N  (10,000  Ibf) 
HMMWV— 22,200  N  (5,000  IbO. 

The  accuracy  associated  with  each  of  these  load 


Time  (s) 

Figure  4.  Shear  stress  versus  time  (shear  annulus). 
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Table  2.  Traction  and  motion  resistance  data  for  the  ivheels! tracks  vehicles. 


Traction 

Tire 

Predicted 

Predicted 

Deiisitx/  Depth  pressure  Measured 

l^ross 

External  resistance 

Date _ (rj  (cm)  (kPa)  (SSM1.0)  %A  (SSM2.0)  %A  Meafiimt  Pmikted  9rA 


a.  5-ton 


28  Jan  88 

-8 

80 

11 

207 

0.36 

0.32 

11 

0.33 

-8 

0.037 

0.033 

-11 

11  Feb  88 

-11 

220 

25 

207 

0.31 

0.37 

20 

0.33 

6 

0.087 

0.082 

-6 

11  Feb  88 

-8 

200 

18 

207 

0.36 

0.36 

0 

0.33 

8 

0.047 

0.061 

30 

10* 

2* 

4* 

8+ 

7+ 

22+ 

b.HMMWV 


19  Jan  88 

-2 

550 

5 

138/152 

0.40 

0.33 

-17 

0.37 

-7 

0.027 

0.00 

-100 

21  Jan  88 

-5 

60 

2-5 

138/152 

0.40 

0.37 

-7 

0.37 

-7 

0.027 

0.019 

-30 

23  Jan  88 

-5 

75 

12 

138/152 

0.42 

0.42 

0 

0.37 

-11 

0.047 

0.070 

49 

26  Jan  88 

-5 

— 

9 

138/152 

0.29 

— 

— 

0.37 

27 

0.067 

— 

— 

28  Jan  88 

-8 

80 

11 

138/152 

0.32 

0.42 

31 

0.37 

16 

0.027 

0.065 

30 

10  Feb  88 

-8 

220 

31 

138/152 

0.29 

0.55 

90 

0.37 

28 

0.117 

0.176 

50 

10  Feb  88 

-7 

250 

22 

138/152 

0.38 

0.53 

40 

0.37 

-2 

0.077 

0.116 

51 

13  Feb  88 

-10 

290 

30 

138/152 

0.35 

0.57 

63 

0.37 

6 

0.127 

0.140 

10 

26  Jan  88 

-5 

148 

18 

138/152 

0.42 

0.47 

12 

0.37 

-12 

0.047 

0.111 

136 

13  Feb  88 

-14 

240 

17 

138/152 

0.28 

0.51 

82 

0.37 

33 

0.097 

0.092 

-5 

33*  5*  33* 

37t  15+  47+ 


e.LAV 

12.5R20 


19  Jan  88 

-2 

550 

5 

207 

0.28 

0.27 

4 

0.35 

25 

0.025 

0.00 

-100 

103 

0.35 

0.41 

17 

0.37 

6 

0.021 

0.00 

-100 

21  Jan  88 

-5 

60 

2-5 

207 

0.45 

0.30 

-33 

0.35 

-22 

0.035 

0.007 

-80 

103 

0.45 

0.41 

-7 

0.37 

-18 

0.021 

0.004 

-81 

23  Jan  88 

-5 

75 

12 

207 

0.34 

0.31 

-9 

0.35 

3 

0.035 

0.026 

-26 

103 

0.45 

0.42 

-7 

0.37 

-18 

0.021 

0.014 

-33 

28  Jan  88 

-8 

80 

11 

207 

0.43 

0.31 

-8 

0.35 

-19 

0.045 

0.024 

-47 

103 

0.48 

0.42 

12 

0.37 

-23 

0.051 

0.013 

-75 

28  Jan  88 

-8 

80 

11 

207 

0.31 

031 

0 

0.35 

13 

0.045 

0.024 

-47 

103 

0.46 

0.42 

-9 

0.37 

-20 

0.051 

0.013 

-75 

28  Jan  88 

-8 

150 

11 

207 

0.28 

0.32 

14 

0.35 

27 

0.045 

0.028 

-38 

103 

0.31 

0.44 

39 

0.37 

19 

0.071 

0.015 

-79 

13  Feb  88 

-10 

290 

30 

207 

0.30 

0.35 

18 

0.35 

18 

0.115 

0.060 

-48 

103 

0.36 

0.49 

36 

0.37 

3 

0.051 

0.032 

-37 

13  Feb  88 

-12 

250 

23 

207 

0.25 

0.34 

38 

0.35 

43 

0.095 

0.052 

-45 

103 

0.25 

0.47 

88 

0.37 

48 

0.081 

0.028 

-65 

11.0R16 

13  Feb  88 

-10 

290 

30 

289 

0.33 

0.31 

6 

0.32 

3 

0.116 

0.085 

-27 

165 

0.40 

0.38 

-5 

0.35 

-12 

0.059 

0.047 

-20 

13  Feb  88 

-12 

250 

23 

289 

0.25 

0.30 

20 

0.32 

28 

0.096 

0.071 

-26 

165 

0.36 

0.37 

3 

0.35 

3 

0.059 

0.041 

-30 

15  Feb  88 

-5 

— 

— 

289 

0.42 

— 

— 

— 

— 

0.016 

— 

— 

165 

0.45 

— 

— 

— 

— 

0.019 

— 

— 

15  Feb  88 

-5 

270 

30 

289 

0.31 

0.30 

3 

0.32 

7 

0.076 

0.089 

17 

165 

0.38 

0.38 

0 

0.35 

-8 

0.069 

0.051 

-26 

15  Feb  88 

-5 

270 

30 

289 

0.30 

0.30 

0 

0.32 

-7 

0.026 

0.089 

242** 

165 

0.43 

0.38 

-12 

0.35 

-8 

0.029 

0.051 

76 

8*  4*  -M* 

24+  20+  38+ 
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Table  2  (cont'd).  Traction  and  motion  resistance  data  for  the  ivheels! tracks  vehicles. 


Traction 


Dale 

Temp. 

(•a 

Density 

(kg/m') 

Depth 

kill) 

Tire 

pressure 

(kPa) 

Measured 

gross 

Predicted 

gross 

(SSM1.0) 

%A 

Predicted 

gross 

(SSM2.0) 

%A 

External  resistance 
Measured  Predicted 

%A 

d.HEMTT 

21  Jan  88 

-5 

60 

2-5 

241/275 

0.32 

0.29 

-10 

0.33 

3 

0.018 

0.006 

-67 

138/207 

0.31 

0.33 

10 

0.33 

-6 

0.014 

0.005 

-64 

23  Jan  88 

-2 

75 

10 

241/275 

0.31 

0.29 

-6 

0.33 

6 

0.008 

0.019 

137 

138/207 

0.33 

0.35 

6 

0.33 

0 

0.024 

0.015 

-38 

28  Jan  88 

-8 

80 

11 

241/275 

0.35 

0.29 

-17 

0.33 

-6 

0.028 

0.021 

-25 

138/207 

0.37 

035 

-6 

0.33 

-11 

0.034 

0.017 

-50 

11  Feb  88 

-11 

220 

25 

241/275 

0.36 

0.31 

-13 

0.33 

8 

0.028 

0.047 

67 

138/207 

0.40 

0.38 

-5 

0.33 

-17 

0.054 

0.032 

^1 

11  Feb  88 

-8 

200 

18 

241/275 

0.26 

0.30 

16 

0.33 

27 

0.028 

0.035 

25 

138/207 

0.37 

0.37 

0 

0.33 

-11 

0.054 

0.024 

-56 

-3’ 

-1’ 

11’ 

10+ 

12+ 

67+ 

e. 

susv. 

21  Jan  88 

-5 

60 

2-5 

_ 

0.51 

2.24 

339 

0.54 

6 

0.060 

^.002 

-97 

26 Jan  88 

-5 

148 

9 

— 

0.57 

228 

300 

0.54 

-5 

0.090 

0.005 

-94 

28  Jan  88 

-10 

150 

11 

— 

0.59 

2.28 

286 

0.54 

8 

0.100 

0.006 

-94 

10  Feb  88 

-8 

220 

31 

— 

0.63 

2.36 

275 

0.54 

-14 

0.100 

0.017 

-83 

10  Feb  88 

-7 

250 

22 

__ 

0.63 

2.34 

271 

0.54 

-14 

0.100 

0.011 

-89 

13  Feb  88 

-10 

290 

30 

— 

0.63 

2.37 

276 

0.54 

-14 

0.080 

0.014 

-83 

13  Feb  88 

-12 

250 

23 

— 

0.62 

2.34 

277 

0.54 

-13 

0.070 

0.012 

-83 

289’ 

-9’ 

-89’ 

22+ 

7+ 

6+ 

f.  M113 

22  Jan  88 

-5 

160 

7 

_ 

0.36 

0.72 

100 

0.42 

17 

0.014 

0.006 

57 

23 Jan  88 

-5 

75 

12 

_ 

0.30 

0.71 

140 

0.42 

40 

0.004 

0.008 

100 

28  Jan  88 

-8 

80 

11 

— 

0.37 

0.71 

92 

0.42 

8 

0.034 

0.008 

-76 

10  Feb  88 

-8 

220 

31 

— 

0.37 

0.75 

103 

0.42 

14 

0.084 

0.024 

-71 

10  Feb  88 

-7 

250 

22 

— 

0.30 

0.75 

150 

0.42 

4 

0.044 

0.016 

-64 

13  Feb  88 

-10 

290 

30 

— 

— 

0.76 

— 

0.42 

— 

— 

0.019 

_ 

13  Feb  88 

-12 

250 

23 

— 

0.41 

0.75 

83 

0.42 

2 

0.084 

0.017 

-80 

111’ 

14’ 

25+ 

13+ 

70+ 

g.  Bradley 

21  Jan  88 

-5 

60 

2-5 

— 

0.36 

0.67 

86 

0.42 

17 

0.011 

0.001 

-90 

23  Jan  88 

-2 

75 

10 

— 

0.31 

0.68 

119 

0.42 

35 

0.011 

0.005 

-54 

28  Jan  88 

-8 

80 

11 

— 

0.35 

0.68 

95 

0.42 

20 

0.031 

0.005 

-84 

12  Feb  88 

-16 

230 

22 

— 

0.45 

0.70 

56 

0.42 

-7 

0.031 

0.011 

-92 

12  Feb  88 

-12 

260 

16 

— 

0.36 

0.70 

94 

0.42 

17 

0.061 

0.008 

-87 

90*  16’ 

20t  13+ 


-82’ 

15+ 


’  Average. 

+  Standard  deviation. 

”  Not  included  in  average  and  standard  deviation  calculation. 


Figure  5.  Shear  strength  versus  normal 
stress  (shear  annulus). 


Table  3  Shear  annulus  test  results  (r  values  [kPa] 
listed). 

Average 

_ g  (kPa) _  density 


Date 

Time 

7.03 

14.06 

21.08 

28.11 

35.14 

(kgjm^) 

22  Jan  88 

14:32 

3.31 

5.65 

7.45 

8.27 

9.45 

80 

23  Jan  88 

14:10 

3.79 

5.93 

7.45 

7.93 

8.89 

60 

23  Jan  88 

14:58 

3.52 

4.83 

7.24 

8.69 

9.24 

60 

26  Jan  88 

11:40 

3.72 

4.69 

6.00 

7.03 

8.07 

130 

26  Jan  88 

13:05 

2.28 

3.45 

4.69 

5.72 

7.58 

130 

26  Jan  88 

13:13 

3.17 

4.34 

5.79 

6.48 

7.86 

130 

28  Jan  88 

10:45 

2.83 

3.79 

5.58 

7.10 

8.76 

150 

28  Jan  88 

10:59 

3.72 

5.65 

6.14 

7.45 

8.34 

150 

n  Feb  88 

15:17 

2.96 

4.00 

5.72 

6.83 

8.48 

200 

13  Feb  88 

9:19 

3.03 

4.83 

5.58 

7.38 

8.27 

220 

13  Feb  88 

9:32 

2.83 

3.86 

4.69 

7.10 

8.27 

220 

13  Feb  88 

10:18 

3.03 

4.27 

5.86 

7.24 

8.41 

230 

13  Feb  88 

10:33 

3.03 

4.34 

5.93 

7.10 

8.62 

230 

13  Feb  88 

12:27 

3.17 

4.41 

5.86 

7.17 

8.69 

270 

13  Feb  88 

12:35 

2.90 

4.48 

6.07 

7.24 

8.48 

270 

13  Feb  88 

13:04 

3.17 

4.27 

5.93 

7.45 

8.69 

239 

15  Feb  88 

14:32 

2.90 

4.27 

5.58 

7.45 

8.55 

250 

17  Feb  88 

14:09 

2.69 

3.72 

4.69 

6.14 

7.58 

280 

18  Feb  88 

11:34 

3.59 

4.83 

5.31 

6.07 

6.96 

260 

18  Feb  88 

14:53 

3.59 

4.62 

4.76 

5.31 

5.79 

280 

18  Feb  88 

14:58 

3.52 

3.93 

4.69 

5.10 

5.72 

280 

18  Feb  88 

15:20 

3.10 

3.79 

4.62 

5.31 

5.93 

260 

28  Feb  88 

15:27 

3.45 

3.93 

4.76 

5.24 

5.79 

260 

25  Feb  88 

11:22 

3.17 

4.14 

4.76 

6.00 

7.45 

180 

25  Feb  88 

11:30 

2.90 

4.48 

5.65 

6.83 

8.21 

180 

25  Feb  88 

13:37 

3.31 

4.62 

6.14 

7.03 

8.21 

190 

25  Feb  88 

13:45 

3.59 

4.83 

5.93 

6.89 

8.14 

190 

2  Mar  88 

15:32 

3.24 

3.93 

4.07 

4.34 

4.62 

220 

2  Mar  88 

15:38 

3.31 

3.31 

4.00 

4.76 

5.58 

220 

3  Mar  88 

8:47 

2.90 

4.14 

5.52 

7.17 

8.48 

240 

3  Mar  88 

8:53 

2.90 

4.21 

5.10 

6.83 

3.21 

240 

3  Mar  88 

13:54 

359 

4.83 

5.38 

6.83 

8.21 

240 

3  Mar  88 

13:59 

4.27 

5.10 

6.14 

7.17 

8.34 

210 

4  Mar  88 

10:41 

4.14 

4.62 

5.03 

5.86 

6.96 

190 

4  Mar  88 

10:46 

4.07 

455 

5.79 

6.41 

7.38 

190 

16  Mar  88 

15:37 

3.31 

3.93 

4.55 

4.76 

5.38 

200 

16  Mar  88 

15:42 

3.10 

3.86 

5.03 

4.83 

5.17 

200 

17  Mar  88 

9:12 

3.10 

3.65 

4.27 

4.27 

5.17 

220 

17Mar88 

9:18 

3.17 

4.21 

4.62 

5.24 

558 

229 

Total  Avg. 

3.24 

4.34 

5.45 

6.48 

7.52 

cells  is  ±  1  /2%  of  full  scale.  Since  accuracy  or  error 
is  based  on  the  full-scale  or  maximum  capacity, 
the  largest  error  will  be  observed  in  the  smallest 
measurement  made  with  a  particular  load  cell. 
Thus,  measurements  for  the  HMMWV  are  within 
1 00  N  (24  Ibf) .  Accuracy  is  within  200  N  (50  IbO  for 
the  44,500-N  load  ceU  and  within  400  N  (100  Ibf) 
for  the  89,000-N  load  cell.  The  test  vehicle  data 
were  all  reported  to  the  nearest  0.01  W,  a  value 
that  is  considerably  larger  than  the  limits  placed 
on  accuracy  by  the  load  cells. 

The  torque  cell  installed  on  the  shear  annulus 
device  had  a  range  of  22.5  N  m  (200  in.  Ibf).  This 
cell  was  calibrated  before  testing  and  several 
times  throughout  the  test  period  to  ensure  that  it 
had  not  been  damaged.  The  normal  range  of 
measured  values  is  from  0  at  the  low  end  to 
around  1 1 .3  N  m  (100  in.  Ibf)  at  the  high  end.  The 
smallest  value  of  torque  that  can  be  read  by  the 
data  acquisition  system  is  approximately  0.()05  N 
m  (0.04  in.  Ibf).  This  is  greater  resolution  than  we 
can  expect  from  the  torque  cell  itself,  since  the  cell 
was  calibrated  using  a  meter  with  an  accuracy  of 
only  0.11  N  m  (1.0  in.  Ibf).  This  results  in  an 
acciuacy  of  about  1  /2%  of  full  scale. 

Snow  conditions 

For  each  test  with  any  vehicle  or  for  tests  with 
the  shear  annulus  device,  KRC  and  CRREL  both 
made  a  set  of  snow  characterizations.  This  in¬ 
cluded  measuring  density  throughout  the  snow's 
depth,  determining  its  temperature,  crystal  size 
and  type,  noting  the  location  of  discontinuities 
such  as  ice  and  frozen  snow  layers,  and  visually 
observing  its  free  water  content.  The  results  of  the 
measurements  made  by  KRC  are  found  in  Ap¬ 
pendix  D.  Table  4  is  a  condensed  listing  of  the 
same  results. 
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Table  4.  Summary  of  snow  characteristics. 


Date 


Total  A^’g-  moio 

depth  Visual  teiiifieralure  deiisili/ 
Time  Site*  (cm)  uvtuessf  (“T)  (kg/m^) 


21  Jan  88 

16;00 

Slopes 

56.0 

dry 

-1.0 

230 

22  Jan  88 

15;00 

Rink 

14.0 

dry 

-5.0 

80 

23  Jan  88 

14;25 

Rink 

10.0 

dry 

-1.0 

<60 

26  Jan  88 

14;00 

Runway 

11.0 

dry 

-10.5 

130 

28  Jan  88 

11;20 

Rink 

11.0 

dry 

-8.0 

150 

10  Feb  88 

12;00 

Rinki 

17.5 

dry 

-7.0 

200 

10  Feb  88 

12:00 

Rink2 

27.0 

dry 

-5.0 

210 

11  Feb  88 

16:00 

Rink2 

31.0 

dry 

-8.3 

200 

11  Feb  88 

16:15 

Rinki 

17.5 

dry 

-7.7 

190 

12  Feb  88 

13:10 

Rink2 

29.0 

dry 

-12.3 

210 

12  Feb  88 

13:20 

Rinki 

16.0 

dry 

-12.5 

220 

13  Feb  88 

9:50 

Rinki 

17.0 

dry 

-12.0 

220 

13  Feb  88 

10:51 

Rink2 

30.0 

dry 

-9.0 

230 

13  Feb  88 

11:00 

Slope's 

83.0 

dry 

— 

320 

13  Feb  88 

11:18 

Texasl 

39.0 

dry 

-7.0 

230 

13  Feb  88 

11:18 

Texas 

76.0 

dry 

-7.0 

260 

13  Feb  88 

12:51 

Sorun 

32.0 

dry 

-8.0 

270 

13  Feb  88 

13:17 

Norun 

23.0 

dry 

-10.3 

230 

15  Feb  88 

14:44 

Rink2 

31.0 

dry 

-5.0 

250 

16  Feb  88 

10:00 

Slope's 

41.0 

dry 

-3.2 

250 

17  Feb  88 

14:30 

Rinki 

17.0 

dry 

-3.7 

280 

18  Feb  88 

11:40 

Rink2 

28.0 

dry 

-2.0 

260 

18  Feb  88 

15:09 

Rinki 

12.5 

wet 

-1.5 

280 

18  Feb  88 

15:33 

Rink2 

20.0 

wet 

-1.3 

260 

25  Feb  88 

11:45 

W/T  Field 

19.0 

dry 

-5.0 

180 

25  Feb  88 

13:55 

Road 

30.0 

dry 

-6.0 

190 

2  Mar  88 

16:00 

W/T  Field 

10.0 

wet 

0.0 

220 

3  Mar  88 

9:05 

W/T  Field 

11.0 

dry 

-9.5 

240 

3  Mar  88 

14:10 

W/T  Field 

12.0 

wet 

-1.0 

210 

4  Mar  88 

11:00 

W/T  Field 

7.0 

dry 

-2.5 

190 

16  Mar  88 

15:45 

W/T  Field 

21.0 

wet 

-0.2 

200 

17  Mar  88 

9:35 

W/T  Field 

25.0 

dry 

-3.3 

220 

’Slopes;  Scriesofhillsofvariablegradeconstructed for vehieleslopcclimbing evalu¬ 
ation. 

Rink;  Large  level  (graded)  area  of  soil,  devoid  of  vegetation,  often  fltxxled  to 
produce  ice  sheet  in  winter.  Rinki  designates  theeast  half  and  Rink2  the  west 
half  of  the  whole  graded  area. 

Runway;  SouthwestendofNE-SW  runway  at  Houghton  County  Airport.  Norun  des¬ 
ignates  the  north  portion  and  sorun  the  south  portion  of  this  section  of  the 


runway. 

Texas;  A  large  field  with  short,  naturally  occurring  grass  cover — surface  is  un¬ 
graded  but  quite  level.  Texasl  designates  a  portion  of  the  Texas  field  where 
vegetation  was  stripped  from  the  surface. 

W/T  Field;  Small  field  cleared  of  old  snow  located  near  KRC  buildings. 

Road;  Access  road  leading  to  W/T  Field, 
t  Determined  by  visual  observation. 


Each  set  of  data  in  Appendix  D  contains  the 
depth  from  the  snow  surface  at  which  a  measure¬ 
ment  was  made.  Densities  are  reported  at  depth 
increments  of  about  0.03  m  because  of  the  size  of  the 
sampling  device.  Ice  and  frozen  snow  layers  are 
indicated  when  present. 

As  can  be  seen  from  the  data,  a  range  of  densities 
was  tested  over  the  period.  Most  of  the  high  density 


snows  resulted  from  their  being  wind  blown;  their 
crystal  size  was  relatively  small. 

ANALYSIS 

The  main  emphasis  of  this  study  was  to  validate 
SSMl  .0  in  its  entirety.  Most  of  our  time  and  effort 
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during  the  first  field  season,  however,  centered  on 
the  traction  algorithm  and  evaluation  of  different 
means  of  obtaining  the  snow  strength  parameters 
c  and  (j). 

Determination  of  snow  strength  parameters 

The  SSM 1 .0  uses  eq  A4  to  generate  snow  strength 
properties  for  input  to  the  traction  calculation  in  eq 
A14.  As  stated  previously,  eq  A4  stems  from  back- 
calculation  of  c  and  (|)  from  actual  mobility  data 
collected  with  the  CIV.  Since  these  c  and  ())  vs 
relationships  resulted  from  averaging  very  scat¬ 
tered  data,  we  were  unsure  of  their  accuracy  and 
perhaps  even  their  validity.  Within  this  study  we 
had  several  means  of  comparing  calculated  values 
of  c  and  4)  with  eq  A4,  including  the  results  of  shear 
annulus  tests,  vehicle  traction  measured  at  two  tire 
deflections  with  some  of  the  test  vehicles,  and  ad¬ 
ditional  traction  measurements  at  two  tire  deflec¬ 
tions  with  the  CIV. 


The  shear  annulus  produced  c  and  ^  values  as 
described  above.  This  results  in  a  c  and  ())  value  for 
every  test  set  (multiple  normal  loads)  (Table  3). 
Plots  of  these  data  against  showed  significant 
scatter,  as  was  also  the  case  for  the  CIV  data  used  to 
produce  eq  A4.  An  alternative  method  of  using  the 
shear  annulus  results  to  obtain  snow  strength  val¬ 
ues  was  then  attempted  to  see  if  clearer  trends 
appeared.  In  this  approach,  we  plotted  measured 
shear  stress  against  p^  for  each  level  of  normal 
load  (Fig.  6).  A  least  squares  regression  routine  was 
applied  to  these  data,  resulting  in  a  series  of  equa¬ 
tions  for  T^as  a  function  of  p^.  TTiese  equations  were 
then  used  to  produce  a  shear  stress  vs  normal  stress 
plot  for  various  densities  (Fig.  7).  Values  of  c  and  4), 
as  a  function  of  p^,  can  then  be  obtained  from  Fig¬ 
ure  7.  This  has  the  effect  of  reducing  the  scatter  in 
the  results  by  recognizing  that  the  field  data  proba¬ 
bly  contain  a  degree  of  unidentifiable  variability. 
Rather  than  to  carry  this  variability  on  into  the  cal- 


(7.  a  =  7.03  kPa. 


b.a  =  14.07  kPa. 


Figure  6.  Shear  stress  versus  initial  snow  density  for  various  normal  loads  (shear  annulus). 
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Shear  Stress  (kPa)  Shear  Stress  (kPa)  Shear  Stress  (kPa) 


Figure  8.  Relationship  between  snoiv  strength  and  initial  density  front  the  shear  annulus  device. 


culations,  this  scheme  removes  the  scatter  at  the 
outset. 

The  results  of  this  exercise  are  shown  in  Figure 
8  and  can  be  given  mathematically  by  the  expres¬ 
sions 

())  =  12.1  -0.0156  p„ 

(5) 

c  =  2.14. 

where  c  is  in  kilopascals,  is  in  kilograms  per 
cubic  meter  and  (j)  is  in  degrees. 

We  indicated  earlier  that,  using  eq  A3,  c  and  ^ 
could  be  calculated  directly  from  data  obtained 
with  a  vehicle.  This  is  accomplished  by  obtaining 
two  values  of  T  through  a  variation  of  either  A  or 
W or  both.  (It  is  usually  more  expedient  to  change 
A  by  altering  tire  pressure,  while  adding  and  sub¬ 
tracting  weight  is  cumbersome.)  The  two  equa¬ 
tions  so  produced  can  then  be  solved  simultane¬ 


ously  for  c  and  ()).  The  CIV  was  operated  at  the  tire 
pressures  of  179  and  103  kPa  (26  and  15  Ib/in.^) 
during  this  test  series  to  allow  determination  of 
snow  strength  from  vehicle  test  results. 

As  with  the  shear  annulus,  two  methods  of  re¬ 
lating  these  data  to  p^  are  possible.  The  first  is  to 
relate  directly  the  individual  values  of  c  and  (()  to 
initial  density,  just  as  was  done  originally  to  pro¬ 
duce  eq  A4.  The  alternate  approach,  as  used  above 
for  the  shear  annulus  device  and  which  will  be 
followed  here  as  well,  attempts  to  remove  experi¬ 
mental  variability  first.  By  determining  linear  best- 
fit  relationships  for  the  vehicle-produced  gross 
tractive  stress  x  vs  p  data  for  each  normal  stress 

g 

level  (Fig.  9),  arid  then  plotting  x  against  normal 
stress  by  solving  these  relationships  for  specific 
values  of  p^^(Fig.  10),  we  obtain  c  and  <|»  as  a  function 
of  density  (Fig.  11).  These  relationships  are  given 
by 
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a.  a  =  179  kPa. 


b.  a  =  103  kPa. 

Figure  9.  Shear  stress  ivrsus  initial  stunv  densitif  for  various  norma!  loads  ( tire). 


^  (degrees)  c  (kPa)  tg ,  Shear  Stress  (kPa) 


c  =  40.67  +  0.0163po 
<t)  =  13.15+ 0.0035  Po 


(6) 
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Comparison  of  c  and  ([)  given  by  eq  A4  (f  rom  Cl  V 
data  obtained  before  1986  [SSMl.O]),  eq  5  (from 
shear  annulus  data)  and  eq  6  (from  1988  CIV  data) 
is  made  in  Figure  1 2.  It  can  clearly  be  seen  that  there 
is  little  agreement  between  these  expressions. 

One  possible  reason  for  the  difference  between 
the  CrV  and  shear  annulus  relationships  is  the 
difference  in  magnitude  of  the  normal  contact 
pressures;  the  CIV  at  the  two  tire  pressures  tested 
had  an  average  contact  pressure  of  184.4  and  201 .6 
kPa,  while  the  shear  annulus  varied  between  7  and 
35.1  kPa.  Along  these  same  lines,  the  shear  rate  for 
the  shear  annulus  is  0.083  m/s,  while  the  tire  oper¬ 
ates  with  shear  rates  of  0  to  greater  than  5  m  /  s,  with 
the  peak  traction  value  occurring  between  0.45  and 
2.2  m/s.  As  noted  by  others  (e.g.,  de  Montmollin 
1982,  McClung  1977),  snow  is  a  very  rate  depend¬ 
ent  material.  The  most  likely  explanation,  how¬ 
ever,  centers  on  the  shear  surfaces;  the  annulus  had 
a  smooth  rubber  face,  while  the  tire  had  a  distinct 
tread  pattern.  This  allows  the  tire  to  utilize  the 
shear  strength  of  the  snow,  while  the  annulus  prin¬ 
cipally  only  engaged  the  surface  friction  character¬ 
istics  of  the  snow.  This  is  very  apparent  in  Figure 
12,  where  the  cohesion  term  is  a  constant  and  the  ()> 
term,  really  a  coefficient  of  friction  v  (v  =  tanc))), 
shows  a  strong  variation  with  initial  snow  density 
for  the  annulus.  It  was  also  very  apparent  in  the 
field  that  the  annulus  engages  only  the  surface 
properties  of  the  snow,  whereas  the  tire  engaged  a 
finite  thickness  of  the  snow  in  the  process  of  shear¬ 
ing.  Because  of  these  gross  differences  in  the  me¬ 
chanics  of  loading  the  snow,  one  should  not  expect 
too  many  similarities  between  the  annulus  and  a 
tire. 

Although  theSSMl  .0  relationships  (eq  A4)  were 
generated  with  CIV  data,  these  c  and  ^  values  were 
established  per  test,  rather  that  via  the  alternate 
method  described  above.  In  part,  this  may  explain 
why  the  new  CIV  expression  (eq  6)  and  the  SSMl  .0 
expressions  (eq  A4)  are  dissimilar. 

Another  possible  conclusion  that  could  be 
reached  from  Figure  12  is  that  c  and  $  as  measured 
here  with  a  tire  are  not  really  functions  of  initial 
density,  or  that  some  other  basic  snow  properties 
exert  a  much  stronger  influence. 

In  studying  the  mechanism  of  a  tire  producing 
traction  on  snow,  and  thinking  about  our  attempt 
to  relate  c  and  ^  derived  from  tire  traction  values  to 
initial  snow  density,  it  soon  became  obvious  to  us 


why  we  were  not  seeing  clear  trends.  The  bulk  of 
the  tractive  force  developed  by  a  tire  or  track  is 
located  in  the  "flat"  portion  of  the  contact  patch. 
The  leading  portion  of  the  snow/tire  or  snow/ 
track  interface  is  usually  smaller  than  the  flat  area, 
and,  because  of  its  curved  geometry,  is  aligned 
more  to  compact  the  snow  than  to  provide  forward 
thrust.  In  the  flat,  horizontal  portion  of  the  contact 
patch,  the  snow  has  been  fully  compacted.  Its  prop¬ 
erties,  especially  its  shear  strength,  are  vastly  dif¬ 
ferent  from  initial  conditions.  Additionally,  since 
we  originally  were  attempting  to  relate  c  and  (])  to 
density,  we  noted  that  the  snow  density  resulting 
from  vehicle  compaction  was  close  to  a  constant  for 
most  vehicles  operating  in  shallow  snow,  varying 
between  400  and  500  kg/m^. 
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These  revelations  suggested  that  it  would  be 
wise  to  either  search  for  other  basic  initial  snow 
properties  that  influence  the  compacted  snow's  c 
and  (]),  or  to  attempt  to  establish  constants  for  the 
values  of  compacted  snow  c  and  (j)  to  replace  ex¬ 
pressions  like  those  in  eq  A4,  5  and  6.  We  hoped 
that  analysis  of  the  traction  data  for  the  wheels/ 
tracks  vehicles  would  give  us  a  clue  as  to  which 
approach  was  most  viable. 

Traction  analysis 

The  simplest  and  probably  most  common  means 
of  reporting  traction  for  a  vehicle  is  by  use  of  a 
coefficient.  Customarily,  this  tractive  coefficient  p 
is  calculated  from  net  traction  measurements.  For 
this  study  we  will  use  gross  traction  and  define  p 
for  a  vehicle  by 

N 

S  % 

|i  -  (7) 

N 

I  ^Vi 
i  =  1 


where  the  summation  is  over  the  driven  tires  or 
tracks  only.  The  measured  gross  traction  columns 
in  Table  2  list  the  average  values  of  p  calculated  for 
all  of  the  test  repetitions  for  each  snow  surface. 

Use  of  p  is  fine  for  comparing  surfaces  and  tires' 
ability  to  support  or  generate  traction,  and  will 
serve  to  allow  comparison  of  the  output  of  SSMl  .0 
with  measured  traction;  however,  by  itself,  it  does 
not  tell  us  much  about  the  mechanics  involved  in 
snow  traction.  Specifically,  we  would  like  to  ex¬ 
plore  further  the  accuracy  of  treating  snow  as  a 
Mohr-Coulomb  material  and  the  relationship  be¬ 
tween  c  and  (})  calculated  from  vehicle  tests  and 
initial  snow  density  p^. 

To  follow  up  on  the  earlier  proposal  that  p^^has 
an  insignificant  effect  on  T.  we  plotted  all  of  the 
traction  data  for  the  wheelsftracks  vehicles,  and  for 
theCrV,  against  p^^for  all  of  the  measured  p  values. 
Figure  13  shows  that  for  all  of  the  vehicles  except 
the  SUSV  and  the  Bradley,  one  could  easily  argue 
that  no  clear  dependence  on  p^  exists.  This  sup¬ 
ports  the  supposition  made  from  examination  of 
the  snow  strength  parameters  data. 


plkg/m^) 


Figure  13.  Tractive  coefficient  versus  initial  snow  density  for  each  vehicle  tested. 
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Figure  14.  Shear  stress  versus  normal 
stress  for  the  23  fan  88  vehicle  tests. 


Assuming  that  the  snow  properties  were  essen¬ 
tially  constant  over  the  span  of  a  single  day,  it 
would  follow  that  c  and  4)  should  remain  fixed 
throughout  the  day  also.  Thus,  if  determination  of 
c  and  <[)  from  vehicle  tests  is  a  valid  means  of 
obtaining  true  material  properties,  data  from  each 
of  the  vehicles  tested  during  a  single  day  should 
give  rise  to  identical  values  of  c  and  <^.  On  23 
January,  gross  traction  data  exist  for  six  vehicles, 
including  three  at  two  contact  pressures.  These 
were  plotted  (Fig.  14)  in  terms  of  normal  (a)  and 
shear  (t  )  stresses.  A  linear  regression  performed 
on  these  data  resulted  in  an  intercept  (c)  of  4.5  kPa 
and  a  slope  (<[»)  of  17.5°  (Fig.  14). 

Recalling  that  SSMl.O  used  only  CIV  data  to 
obtain  c  and  ())  values,  we  also  passed  a  best-fit  line 
through  theCrV  data  by  itself  (Fig.  14).  Clearly,  this 
hne  has  a  lower  slope  (14.5°)  and  a  much  higher 
intercept  (16.1  kPa).  The  significance  of  this  result 
will  be  discussed  later. 

It  is  readily  apparent  in  Figure  14  that  the 
wheeled  vehicle  data  fit  a  straight  line  well,  while 
the  two  data  points  corresponding  to  the  tracked 
vehicles  do  not  lie  along  this  same  line.  This  sug¬ 
gests  that  perhaps  a  nonlinear  expression  describes 
the  relationship  between  normal  load  and  shear 
failure  of  snow. 

Mohr's  original  shear  failure  criterion  for  soils 
stated  only  that  the  shear  stress  at  failure  was  a 
function  of  the  normal  stress  7  on  that  plane. 
(Coulomb  chose  to  define  the  fu:  iction  as  a  straight 
line  for  limited  normal  stress  ranges.)  Since  we  are 
considering  a  wide  range  of  normal  stresses  (val¬ 
ues  of  15  kPa  for  the  SUSV  up  to  265  kPa  for  the 
LAV)  it  is  perfectly  justifiable  (and  still  within 


Mohr-Coulomb  theory  confines)  to  choose  a  non¬ 
linear  function  that  better  fits  these  data. 

Other  single-day  data  were  plotted  in  a  fashion 
similar  to  Figure  14.  Each  resulted  in  very  similar 
distributions  in  the  data.  This  observation,  coupled 
with  our  earlier  proposition  that  the  compacted 
snow  provides  the  bulk  of  the  tractive  force  and 
that  the  compacted  snow  is  not  only  close  to  con¬ 
stant  from  vehicle  to  vehicle,  but  from  day  to  day, 
led  us  to  place  all  of  the  winter's  data  on  one  plot 
(Fig.  15).  A  distinct  trend  is  obvious,  and  one  must 
assume  that  gross  traction  developed  by  tires  or 
tracks  is  principally  a  function  of  normal  stress  on 
the  driven  elements,  regardless  of  the  initial  snow 
properties. 

A  straight  line  placed  through  the  data  in  Figure 
15  yields  c  and  (j)  values  of  4.6  kPa  and  17.9°,  re¬ 
spectively,  with  a  correlation  coefficient  of  0.92  for 
the  fit.  This  is  certainly  a  strong  correlation;  how¬ 
ever,  the  large  number  of  wheeled  vehicle  data  and 
their  distribution  on  the  normal  stress  axis  resulted 
in  a  poor  fit  for  the  tracked  vehicles.  As  noted 
earlier,  there  is  no  reason  not  to  believe  that  a  non¬ 
linear  relationship  governs  shear  failure  in  snow.  A 
power  function  was  thus  chosen  to  model  these 
data.  This  yielded  the  expression 

Tg  =  0.851  o0-823  (8) 

where  o  is  in  kilopascals.  The  correlation  coeffi¬ 
cient  for  eq  8  is  0.97. 

The  scatter  present  in  the  data  displayed  in 
Figure  15,  although  small  by  mobility  test  stan¬ 
dards,  can  perhaps  be  explained  by  several  factors. 
First,  we  have  only  associated  a  single  value  of  a 
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Figure  15.  Shear  stress  versus 
normal  stress  using  all  data. 


with  each  vehicle.  This  value  is  determined  stati¬ 
cally  and  is  an  average  for  the  whole  vehicle,  being 
gross  vehicle  weight  divided  by  the  number  of 
contact  points  and  the  average  contact  area  of  these 
points  as  measured  on  a  hard  paved  surface.  Ob¬ 
viously,  during  testing  on  deformable  surfaces 
under  dynamic  conditions,  both  contact  area  and 
normal  load  may  vary  considerably.  If  an  average 
A  and  a,  for  each  tire  or  track,  were  available  to  as¬ 
sociate  with  the  measured  T  ,  for  each  test,  a  more 

n 

accurate  plot  could  be  produced. 

Additionally,  Figure  1 5  plots  gross  tractive  stress, 
which  is  the  result  of  adding  the  measured  net 
traction  T  with  the  measured  motion  resistance  R  . 

n  s 

Recall  that  these  quantities  are  collected  from  two 
separate  tests,  one  with  the  tractive  elements  driv¬ 
ing,  the  other  with  them  free  rolling.  During  the 
motion  resistance  test,  only  sinkage  attributable  to 
the  vehicle's  normal  stress  occurs,  while  the  trac¬ 
tion  test  may  generate  both  normal  load  sinkage 
and  sUp  sinkage.  Thus,  the  value  of  we  add  to 
to  obtain  T  does  not  account  for  the  added  motion 

g 

resistance  caused  by  slip  sinkage. 

Another  possible  explanation  for  the  variation 
in  measured  values  in  Figure  15  is  that  a  given  test 
may  not  have  been  conducted  in  a  manner  that 
resulted  in  engaging  the  full  available  shear  strength 
of  the  snow.  This  happens  when  a  tire  or  track 
experiences  little  or  no  slippage  (i.e.,  power  limita¬ 
tion)  or  if  only  a  few  tires  on  a  vehicle  slipped  (i.e., 
action  of  differentials).  All  values  of  shear  stress 
below  the  shear  strength  of  snow  are  acceptable. 


Lastly,  as  with  other  geologic  materials,  snow 
properties  vary  in  time  and  space.  The  strength, 
depth,  structure,  water  content,  temperature,  and 
even  type  and  strength  of  supporting  material, 
undoubtedly  were  not  constant  along  a  single  test 
lane  or  from  lane  to  lane.  This  gives  rise  to  vari¬ 
ations  in  sinkage  and  to  the  available  shear  strength. 
Considering  these  factors,  we  would  expect  some 
variation  in  for  any  given  vehicle. 

Traction  model  predictions 

To  test  the  viability  of  SSMl  .0,  the  model  was 
run  for  all  of  the  snow  conditions  at  which  data 
were  collected  with  the  loheels/ tracks  family  of 
vehicles.  In  this  way,  actual  and  predicted  values  of 
traction  could  be  compared.  Table  2  contains  a 
listing  of  the  measured  and  predicted  values  of 
gross  traction.  Another  column  also  shows  the 
degree  to  which  the  predicted  values  vary  from  the 
measured  values.  The  errors  range  from  a  4%  un¬ 
derprediction  for  the  HEMTT  tc  a  298%  overpredic¬ 
tion  of  traction  for  the  SUSV.  In  general,  the  SSMl  .0 
does  a  good  job  of  predicting  for  the  wheeled  ve¬ 
hicles;  an  average  error  of  15%  exists  among  them. 
The  model,  however,  grossly  overpredicts  for  the 
tracked  vehicles,  with  an  average  offset  of  162%. 

While  the  results  of  the  SSMl.O  prediction  are 
encouraging  for  the  wheeled  vehicles,  quite  the 
opposite  can  be  said  for  tracks.  We  recall  that  the 
model  predicts  traction  based  on  application  of  eq 
1,  using  eq  Al,  A2and  A4  (generated  with  the  CIV). 
Returning  to  Figure  14,  we  see  that  using  only  the 
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CIV  data  to  produce  a  straight-line  expression  to 
describe  the  normal  stress-shear  stress  relation¬ 
ship  results  in  grossly  overshooting  the  tracked 
vehicle  data.  This  line,  however,  is  not  too  far  from 
passing  through  the  majority  of  the  wheeled  ve¬ 
hicle  data. 

Clearly,  theCIV-produced  line  has  a  much  higher 
intercept  (c)  and  a  slightly  lower  slope  (4))  than  the 
line  generated  for  all  of  the  vehicle  data.  Thus,  the 
( W  tan(]))  term  from  ec]  A3  should  ■  'iffer  little  for 
either  wheeled  or  tracked  vehicles,  depending  on 
which  (|)  is  used.  The  two  c  values  are  significantly 
different,  leading  to  considerable  differences  in  the 
(Ac)  term  from  eq  A3.  For  wheeled  vehicles,  this 
term  contributes  little  to  gross  traction;  however, 
for  tracked  vehicles,  their  large  contact  area  makes 
this  term  the  dominant  one. 

Obviously,  extrapolation  of  the  CIV  data  down 
to  the  contact  pressure  range  for  tracked  vehicles 
results  in  huge  over-predictions.  Interestingly,  the 
prediction  error  magnitudes  for  the  whole  family 
of  zuheels/tracks  vehicles  is  in  direct  proportion  to 
the  difference  between  their  ground  pressures  and 
the  CIV'S. 

The  SSMl  .0  uses  initial  snow  density  to  calculate 
traction  for  the  leading  tire  and  various  compacted 
snow  densities  (depending  on  the  contact  pressure 
for  that  tire  or  track)  for  trailing  tires  (eq  A14). 
Figure  15  indicates  that  this  scheme  is  not  neces¬ 
sary.  The  snow  that  is  supporting  tractive  forces  is 
essentially  constant  for  all  vehicles  and  conditions 
within  the  range  of  snows  tested  in  this  study. 

In  light  of  the  prediction  errors  associated  with 
SSMl.O  and  the  new  relationship  portrayed  in 
Figure  15,  we  instituted  a  change  in  the  traction 
prediction  algorithm  contained  in  SSMl.O.  The 
updated  model,  SSM2.0  (Appendix  E,  including 
sample  output;  flow  chart  given  in  Fig.  16),  now 
predicts  traction  based  on  eq  8.  Although  this  may 
seem  to  result  in  an  entirely  empirical  relationship 
specific  to  this  study,  we  use  two  points  to  defend 
the  introduction  of  a  new  traction  algorithm  for 
future  use.  First,  the  data  plotted  in  Figure  15  show 
a  very  good  correlation  (correlation  coefficient  of 
0.97)  to  a  power  curve.  The  new  relationship  is 
based  on  the  same  premise  used  in  the  previous 
traction  model,  and  with  the  most  popular  traction 
prediction  method  referenced  in  the  literature  over 
the  past  30  years  (eq  A3).  Additionally,  the  data  in 
Figure  15  are  for  a  wide  variety  of  vehicles — ve¬ 
hicles  with  gross  differences  in  weight,  traction 
elements,  configurations,  etc. — yet  a  strong  corre¬ 
lation  exists.  This  suggests  to  us  that  eq  8  depicts  a 


universal  relationship  for  snow  traction,  rather 
than  a  one-time,  site-specific  correlation. 

Our  second  defense  rests  on  comparing  three 
external  data  sets  (data  generated  in  other  studies) 
with  eq  8.  The  LAV  was  tested  in  sometimes  deep 
alpine  snow  at  Pickle  Meadows,  California,  several 
years  ago.  These  data  (Willoughby,  in  prep.)  show 
an  average  agreement  with  eq  8  of  15%  with  a 
standard  deviation  of  9%.  Also  available  are  mobil¬ 
ity  data  for  the  prototype  Hard  Mobile  Launcher 
tractor  (HMD  (8x8)  from  Great  Falls,  Montana 
(Blaisdell  et  al.  1987).  The  snow  traversed  by  the 
HML  was  wet  and  sometimes  very  shallow.  Agree¬ 
ment  with  eq  8  for  the  HML  was  also  very  good, 
with  an  average  agreement  within  5%  (standard 
deviation  is  21%).  Mobility  data  for  the  Challenger 
C65  rubber-belted  tractor  are  also  compared  here. 
These  data  (Blaisdell  and  Liston  1 989),  for  two  track 
widths,  display  an  average  agreement  of  26%  and 
standard  deviation  of  3%  with  the  model's  predic¬ 
tions.  All  of  these  external  data  (also  shown  in  Fig. 
15)  agree  with  the  ivheels/tracks  data  quite  well, 
despite  being  from  different  test  locations.  Two  of 
these  data  sets  were  also  gathered  with  different 
vehicles  from  those  used  in  the  loheels/ tracks  study 
(the  Challenger  and  HML),  yet  they,  too,  compli¬ 
ment  the  ivheels/tracks  vehicle  data. 

The  second  generation  shallow  snow  model, 
SSM2.0,  uses  the  gross  traction  relationship  (for  a 
single  driven  tire  or  track) 

Tg.  =0.851  Ai(WiMi)°®“3  (9) 

where  W  is  in  kilonewtons  and  A  is  in  square 
meters.  In  SSM2.0,  no  adjustments  need  to  be  made 
for  different  areas  or  percent  of  new  snow  com¬ 
pacted  or  snow  conditions  such  as  density  and 
depth,  as  was  the  case  in  SSMl.O  (see  eq  A14).  In 
essence,  the  result  is  that  gross  traction  on  snow 
becomes  a  constant  for  a  given  vehicle  ground 
pressure,  and  that  differences  in  for  different 
snows  or  different  vehicles  are  the  result  of  differ¬ 
ent  magnitudes  of  R^.  (This  is  limited  to  shallow  dry 
snows  with  densities  of  less  than  500  kg/ m^.) 

Using  the  SSM2.0,  T^was  also  calculated  for  the 
cases  where  actual  wheels/tracks  data  were  taken. 
Table  2  includes  a  column  with  the  new  model 
prediction  and  its  degree  of  accuracy  is  compared 
with  the  measured  values.  For  the  ivheels/tracks 
vehicles,  differences  between  predicted  and  meas¬ 
ured  range  from  -9  to  1 6% .  The  SSM2.0  still  shows 
a  larger  difference  for  the  tracked  vehicles  than  for 
the  wheeled  vehicles,  but  a  much  more  acceptable 
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Figure  1 7.  Noniialized  motion  resistance  as  a  function  of  a  deptli-tinies-densit}/  snozo  parameter. 
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level  of  difference  exists  between  predicted  and 
measured  gross  traction  in  all  cases. 

Resistance  analysis 

As  described  above  (eq  1),  a  measure  of  vehicle 
mobility  can  be  determined  from  the  difference 
between  the  gross  traction  T  and  external  motion 
resistance  R^.  Although  this  study  concentrated 
primarily  on  the  prediction  of  T  ,  the  accurate 
determination  of  is  equally  important.  Motion 
resistance  data  were  used  to  1)  calculate  T  for 
validation  of  the  traction  algorithm,  and  2)  validate 
the  resistance  algorithm  used  in  the  SSMl.O  (or 
SSM2.0,  since  version  2.0  differs  from  version  1 .0 
only  in  the  traction  algorithm). 

As  with  gross  traction,  resistance  was  presented 
in  coefficient  form  (R^/  W),  and  this  coefficient  will 
be  used  to  compare  measured  and  predicted  resis¬ 
tance.  This,  however,  does  not  provide  a  great  deal 
of  insight  into  the  mechanisms  taking  place  that 
give  rise  to  R^.  Our  hope,  once  again,  was  to  de¬ 
velop  a  simple  relationship  between  vehicle  para¬ 
meters  and  snow  characteristics. 

The  majority  of  motion  resistance  for  the  ve¬ 
hicles  used  in  this  study  was  assumed  to  occur  at 
the  front  axle;  the  bulk  of  the  snow  compaction 
takes  place  here.  The  total  sinkage  after  vehicle 
passage  will  be  governed  by  the  greatest  axle  load. 
We  also  note  that  tire  or  track  width  should  also  be 
an  important  factor.  Two  of  the  simply  measured 
snow  properties  that  affect  resistance  are  snow 
density  and  depth.  With  these  ideas  in  mind,  we 
constructed  Figure  17,  where  the  x-axis  contains  a 
combined  snow  parameter  (initial  density  times 
depth)  and  the  i/-axis  displays  the  vehicle  parame¬ 
ter  {R^/  b).  The  snow  depth  included  in  the 

.r-axis  parameter  is  limited  to  compactible  snow 
(i.e.,  p  <  500  kg/m^).  All  of  the  data  are  presented 
except  those  for  the  LAV25  with  11.00R16  tires. 
These  data  did  not  fit  the  general  trend  presented 
by  the  other  vehicles  and  lead  us  to  question  their 
validity.  The  LAV25  with  12.5R20  tires  shows  more 
scatter  than  the  other  vehicle  data  but  follows  the 
same  trend.  Straight  lines  are  fit  to  these  data  for 
each  vehicle  and  the  slopes  and  intercepts  are  pre¬ 
sented  in  Table  5.  In  this  table  it  can  be  seen  that  the 
slope  decreases  as  tire  radius  increases  or  as  gross 
vehicle  weight  increases.  The  intercepts  are  greater 
than  zero  and  are  probably  associated  with  the  tire 
or  track  rolling  resistance  on  compacted  snow  where 
no  sinkage  is  displayed .  Compacted  snow  undoubt¬ 
edly  displays  some  elastic  behavior,  so  some  tires 
and  tracks  probably  generate  a  small  amount  of 
recoverable  deformation  as  they  roll  over  the  snow. 


Table  5.  Results  of  regression  analysis  on 
motion  resistance  data. 


Average 


Vehicle 

Slope 

Intercept 

VJeight  (N) 

tire  radius  (m) 

CIV 

0.00946 

0.0314 

24,696.4 

0.352 

HMMWV 

0.0068 

0.253 

33,450.5 

0.429 

LAV25» 

0.0068 

0.663 

125,483.7 

0.429 

5-ton 

0.0057 

0.253 

105,778.2 

0.554 

Bradley 

0.0039 

0.097 

223,299.6 

0.356‘» 

SUSV 

0.00086 

0.485 

61,340.7 

0.264‘» 

M113 

0.00283 

0.137 

104,087.9 

0.368‘» 

10-ton 

0.0037 

0.331 

268360.1 

0.603 

*  LAV25  with  12.5  R20  tires. 

**  Estimated  radius  at  leading  edge  of  track  (from 
SSMl.O). 


This  rolling  resistance  was  not  measured,  but  is  no 
doubt  greater  than  the  resistance  felt  on  an  unde- 
formable  surface.  The  CIV  can  be  considered  as  a 
one-axle  vehicle  in  this  plot  since  only  the  front 
wheel  resistance  was  measured.  The  hard  surface 
motion  resistance  of  the  front  wheels  of  the  CIV  is 
on  the  order  of  100  N,  which  corresponds  to  a  value 
of  0.031  on  the  y-axis  of  Figure  17;  this  is  almost 
exactly  as  predicted  by  the  regression  analysis 
(0.0314).  These  observations  suggest  that  a  simple 
empirical  algorithm  could  be  developed  with  which 
to  predict  motion  resistance  in  snow. 

Since  a  second  year's  study  is  planned  for  the 
loheelsf tracks  project,  the  resistance  model  was 
developed  no  further  at  this  time.  Instead,  the 
second  year  of  study  will  concentrate  on  under¬ 
standing  motion  resistance  in  snow. 

Resistance  model  predictions 

We  complete  our  analysis  with  a  comparison  of 
the  SSM1.()  resistance  prediction  with  the  resis¬ 
tance  values  measured  in  the  field  tests.  Table  2 
contains  both  measured  and  predicted  (eq  A13) 
resistance,  in  coefficient  form,  and,  as  for  the  trac¬ 
tion  comparison,  the  percent  difference  between 
the  predicted  and  measured  values  was  calculated. 
The  CIV  data  were  not  compared  with  the  model 
since  the  CIV  measures  R^  at  the  front  wheel  only 
and  the  SSMl.O  predicts  whole  vehicle  motion  re¬ 
sistance.  Motion  resistance  is  generally  under-pre- 


23 


dieted  for  the  tracked  vehicles  (by  as  much  as 
nearly  90%  in  one  case),  while  no  clear  trend  is  seen 
for  the  wheeled  vehicles.  The  lack  of  any  clear  trend 
in  wheeled  vehicle  comparisons  and  the  under¬ 
prediction  of  tracked  vehicle  motion  resistance 
suggests  that  the  resistance  algorithm  may  not  be 
appropriate  for  all  vehicles.  The  possibility  of  ex¬ 
cessive  (arbitrarily  set  at  less  than  ±20%)  error  is 
much  greater  for  the  resistance  algorithm  (80% 
chance  of  excessive  error)  than  for  the  new  traction 
algorithm  and  thus  it  is  the  weakest  part  of  SSMl  .0 
orSSM2.0. 


CONCLUSIONS  AND 
RECOMMENDATIONS 

Use  of  the  CIV-generated  snow  strength-initial 
density  relationship  to  predict  gross  traction  (as  is 
done  in  SSMl  .0)  was  found  to  produce  marginal 
correlation  with  measured  values  for  wheeled 
vehicles,  and  dismal  results  for  tracked  vehicles. 
We  explored  the  possibility  of  producing  the  c  and 
(])  vs  density  relationship  with  a  shear  annulus 
device,  and  found  that  it  did  not  even  show  a  trend 
similar  to  that  found  with  the  CIV.  We  concluded 
that  initial  snow  density,  in  the  range  from  50  to  500 
kg/m^,  has  little  or  no  effect  on  the  traction  levels 
generated  by  vehicles. 

By  plotting  traction  for  all  of  the  vehicles,  in 
terms  of  a  failure  shear  stress,  against  normal  pres¬ 
sure  exerted  by  the  vehicles,  we  discovered  a  Mohr- 
type  failure  relationship.  Regression  analysis  re¬ 
sulted  in  a  very  good  fit  for  a  power  function,  and 
predictions  based  on  this  relationship  show  excel¬ 
lent  agreement  with  measured  values  (averaging 
within  7%)  and  hold  considerable  promise  for  the 
future.  Version  2.0  of  the  shallow  snow  model 
contains  this  new  traction  algorithm. 

The  resistance  algorithm  used  in  SSMl.O  and 
SSM2.0  was  compared  with  the  data  collected  and 
poor  agreement  was  observed.  This  is  clearly  the 
weakest  part  of  the  model  and  thus  must  be  im¬ 
proved.  An  initial  analysis  of  the  resistance  data 
was  made  and  it  appears  that  a  simple  empirical  al¬ 
gorithm  may  be  obtainable.  These  data  may  also  be 


used  to  assist  in  the  development  of  a  new  analyti¬ 
cal  model;  both  approaches  should  be  considered. 
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APPENDIX  A:  SHALLOW  SNOW  MOBILITY  MODEL,  VERSION  1.0. 


DESCRIPTION 


Motion  resistance  is  a  function  of  the  parameters  that  affect  vehicle  sinkage  and  the 
tendency  to  displace  snow  laterally.  A  partial  list  includes  the  load,  contact  pressure,  snow 
strength  and  depth,  and  width  of  the  tire  or  track.  During  the  past  30  years,  several  resistance 
models  have  been  proposed  in  the  literature  for  deformable  materials.  With  the  goal  of 
keeping  SSMl.O  simple  (i.e.,  a  model  that  has  a  short  list  of  input  parameters  that  are  easily 
obtained)  and  allowing  it  to  address  a  broad  range  of  vehicle  and  snow  conditions,  these 
resistance  expressions  were  scrutinized.  The  vehicle  data  required  by  each  of  these  expres¬ 
sions  are  similar  from  model  to  model,  and  are  readily  accessible.  The  snow  data  required 
to  process  any  of  these  expressions,  however,  vary  considerably.  The  only  model  that 
requires  snow  data  that  can  be  quickly  and  routinely  obtained  in  the  field  is  that  of  Liston 
(1974). 

Liston  assumes  that  a  hyperbolic  relationship  exists  between  compacting  pressure  and 
volume.  Applying  energetics,  he  then  integrates  between  the  initial  and  final  volumes  of 
snow  involved  in  compaction  to  obtain  the  work  of  snow  compaction.  This  implicitly  as¬ 
sumes  that  there  is  no  bulldozing  (displacement)  of  the  snow.  If  there  is  no  volume  change, 
no  work  is  done.  Finally,  the  work  of  compaction  is  equated  to  motion  resistance  times  the 
horizontal  distance  traveled. 

If  we  assume  that  lateral  flow  of  the  snow  during  compaction  is  insignificant  (i.e., 
compaction  is  confined  to  the  width  of  the  tire  or  track),  volume  change  in  the  snow  can  be 
expressed  in  terms  of  the  sinkage  z.  Further,  if  the  total  mass  of  the  snow  does  not  change 
during  compaction,  then  initial  and  final  volumes  of  snow  can  be  related  to  the  initial  and 
final  densities  of  snow.  We  can  then  write 


Ks  =  Po!  [V(pf  -po)]  ln(pf/po)  - (1/pf) } .  Po  <  Pf  (z  >  0) 
Rs  =  0  p„  =  p,(z  =  0) 


where  p 
b 
h 

Po 

Pf 


z 


tire  inflation  pressure 
maximum  tire  or  track  width 
snow  depth 

initial  snow  density  (prior  to  tire  or  track  passage) 
ending  snow  density  (after  tire  passage) 
sinkage. 


(Al) 


Driving  traction  is  also  a  sum  of  the  interaction  of  many  snow  and  vehicle  parameters. 
Those  that  were  mentioned  above  for  resistance  still  apply,  along  with  more  detailed  features 
of  the  tire  or  track  (e.g.,  tread  pattern,  tire  or  track  "rubber"  compound,  grouser  spacing  and 
height,  grouser  or  tread  geometry).  The  number  of  traction  models  proposed  in  the  literature 
is  fewer  than  is  the  case  for  motion  resistance.  These  models  seem  to  fall  into  two  categories — 
either  they  are  very  simplistic,  lumping  many  parameters  together  into  a  few  constants,  or 
they  are  exceedingly  detailed. 

We  assumed  for  SSMl.O  that  the  gross  tractive  force  T^in  snow  is  developed  by  a  tire  or 
track  through  a  combination  of  frictional  forces  and  engaged  snow  shear  strength  developed 
in  the  tire/snow  interface  region.  The  magnitude  of  the  frictional  force  at  the  interface  is  a 
function  of  the  snow  type  and  the  tire  tread  compound.  Shear  forces  are  developed  by  the 
internal  friction  and  cohesion  of  the  snow. 
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Opting  once  again  to  keep  the  SSMl.O  simple,  our  traction  model  is  based  on  the 
Mohr-Coulomb  relationship  for  soil  shear  failure 

s  =  c  +  a  tancj) ,  (A2) 

where  s  =  shear  strength  at  failure 
a  =  normal  stress 
c  =  cohesion 

<|)  =  internal  angle  of  friction. 

Use  of  this  criterion  for  snow  could  be  questioned,  but  Yong  and  Fukue  (1978)  present  direct 
shear  test  data  for  snow  that  can  be  represented  by  eq  A2  quite  well.  Additionally,  eq  A2  has 
been  used  as  the  basis  for  modeling  traction  for  some  time,  and  the  more  rigorous  models  of 
snow  failure  are  far  too  complicated  for  the  mobility  model  we  sought. 

For  mobility  purposes,  eq  A2  is  rewritten  in  terms  of  forces.  Multiplying  through  by  area, 
we  obtain 


Tg  =  Ac  +  W  tan<t) 


(A3) 


where  T  is  the  highest  possible  gross  traction  attainable  by  a  single  driven  tire  or  track,  or 
the  maximum  shear  strength  available  from  the  snow.  A  is  the  tire  or  track's  contact  area. 
Unfortunately,  the  snow  strength  parameters  c  and  <])  are  not  readily  available  in  all  cases. 
Historically,  c  and  (])  are  determined  from  shear  annulus  (e.g.,  bevameter)  test  results  using 
eq  A2.  It  is  also  possible  to  obtain  c  and  <|)  from  eq  1  and  A3  using  actual  vehicle  tests,  where 

net  traction  and  motion  resistance  data  are  col- 


Po(kg/m3) 


lected  at  two  or  more  different  inflation  pres¬ 
sures  or  at  two  or  more  different  normal  loads,  or 
both. 

Using  eq  A3  for  calculation  of  gross  traction 
naturally  limits  SSMl.O's  predictions  to  those 
cases  where  c  and  (])  data  exist.  Because  of  the 
paucity  of  these  snow  strength  parameters  for 
snow  and  questions  about  the  reliability  of  c  and 
<[)  values  obtained  with  a  shear  annulus,  we 
sought  a  means  of  determining  these  parameters 
from  snow  data  that  were  more  readily  available 
or  predictable. 

Since  initial  snow  density  is  the  most  sophisti¬ 
cated  of  the  simply  and  customarily  measured 
snow  properties,  we  attempted  to  relate  the 
strength  parameters  c  and  (])  to  this  snow  prop¬ 
erty.  Each  set  of  traction  data  available  from  CIV 
experiments  conducted  before  1986  was  used  to 
determine  c  and  (])  values  for  each  traction  test. 
These  values  were  subjected  to  a  least  squares 
regression  analysis  with  initial  snow  density 
as  the  independent  parameter.  Although  con¬ 
siderable  scatter  was  present  in  the  c  and  (() 


Figure  AI.  Relationship  between  snow  versusdensityplots,  the  SSMl.O  used  the  results 
strength  and  initial  density  as  used  in  of  this  regression  analysis,  which  can  be  stated  as 

SSM2 .0.  (and  shown  in  Fig.  Al ) 
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c 


0.66  +  0.0887  po 
17.17 


0  <  Po  <  300 
300  <  Po  <  500 


(A4) 


(t)  =  17.06  -  0.01 63  Po 

where  c  is  in  kilopascals,  p^is  in  kilograms  per  cubic  meter  and  ([)  is  in  degrees. 

The  mobility  expressions  given  by  eq  1,  A1  and  A3  should  be  thought  of  for  a  single  tire 
or  traction  element.  For  a  mobility  model  to  be  flexible,  neither  the  specific  nor  general 
configuration  of  a  vehicle  should  be  Umited  by  the  model.  Being  per  tire  or  per  track,  the 
relationships  for  traction  and  resistance  given  here  are  used  on  the  vehicle's  tires  or  tracks, 
one  at  a  time,  in  the  SSMl  .0  until  all  of  the  traction  elements  have  been  accounted  for.  In  this 
way,  tires  or  tracks  with  different  loads,  inflation  pressures,  sizes,  configurations  (dual  or 
single,  driven  or  free-wheeling),  degree  of  tracking,  and  position  on  the  vehicle  are  all 
accommodated  for  with  one  set  of  equations.  Placed  in  a  loop  in  SSMl  .0,  these  expressions 
are  used  for  each  station,  and  a  running  sum  for  the  whole  vehicle  is  accumulated  to  produce 
a  measure  of  the  net  performance  of  the  vehicle.  A  station  is  defined  here  as  a  transverse 
section  of  the  vehicle  including  a  single  set  of  wheels  (single  or  dual  tire)  or  tracks  on  each 
side  of  the  vehicle  (i.e.,  the  vehicle  is  assumed  to  be  symmetric  about  its  longitudinal 
centerline). 

Clearly,  to  apply  the  traction  and  resistance  equations  above  to  even  a  single  tire  or  track, 
it  is  necessary  that  we  calculate  several  parameters  (as  they  are  not  readily  known  or 
measurable  before  or  after  the  vehicle  is  in  the  field).  In  the  determination  of  T^,  it  is  necessary 
to  know  the  area  of  contact  A  between  the  tire  or  track  and  the  snow .  As  will  be  apparent  later, 
the  sinkage  of  the  tire  or  track  is  required  to  find  A.  To  calculate  R^,  initial  and  final  densities 
associated  with  the  passage  of  a  tire  or  track  are  necessary.  In  the  following  paragraphs,  we 
first  describe  how  these  parameters  are  determined  and  then  proceed  to  show  how  the 
mobility  equations  are  adapted  for  application  to  a  whole  vehicle  in  SSMl  .0. 

Contact  area  is  determined  using  a  simplistic  approach  where  A  is  broken  down  into  two 
regions.  The  flat  area  directly  under  the  tire  or  track,  Aj,  is  calculated  from 


A2  =1  I 

)  bf  (track)  ( 


(A5) 


where  t  is  the  track  length.  The  tire's  sidewall  strength  has  been  ignored.  The  remainder  of 
the  area  is  curved  and  occurs  as  a  result  of  the  tire's  or  track's  sinkage  into  the  snow.  Using 
a  geometric  approach,  this  area  can  be  found  from  the  expression 


Ai 


=  br  arccos 


r-z 


(A6) 


where  r  is  the  tire's  radius  or  a  theoretical  radius  for  a  tracked  vehicle  (Fig.  A2).  The  total 
contact  area  A  is  the  sum  of  these  two  components;  obviously,  tires  or  tracks  that  exhibit  no 
sinkage  have  no  curved  contact  area  component. 

Using  eq  A5  and  A6  to  calculate  A,  we  must  know,  at  each  station,  the  vehicle  parameters 
p,  W  (for  wheeled  vehicles  only),  f  (for  tracked  vehicles  only),  and  b  and  r  of  the  tire  or  track. 
For  the  snow,  it  is  required  that  we  know  z.  To  determine  sinkage,  it  is  assumed  that 
compaction  (within  the  realm  of  loads  that  are  most  common  for  vehicles)  only  occurs  in  the 
vertical  dimension.  First,  we  claim  that  the  maximum  sinkage  that  will  occur  as  the 
result  of  vehicle  passage  can  be  calculated  from 
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Zfinal  -  o  ( 1  / 


where  =  depth  of  the  virgin  snow 

pj  =  maximum  or  final  density  (Yong  and  Fukue  1977)  in  any  rut  following  vehicle 
passage. 

It  is  reasonable  to  believe  that  z^.^^jand  thus  p^will  occur  under  the  tire  or  track  applying  the 
largest  ground  pressure;  p^is  thus  determined  in  the  SSMl.O  based  on  the  highest  pressure 
exerted  by  any  tire  or  track  on  the  vehicle.  Four  major  categories  were  defined  (based 
on  experience)  as  follows: 


Max.  ground  pressure  Final  density 

(kPa)  (kg/m^) 

<  210  500 

211-350  550 

351-700  600 

>  701  650 

Intermediate  values  of  sinkage  z,  for  tires  or  tracks  with  inflation  or  contact  pressures  (p) 
less  than  were  determined  by  applying  the  ratio  of  to  Since  the 

pressure-sinkage  relationship  is  obviously  not  linear  for  compressible  snow  (p^  <  500  kg/ 
m^)  (Abele  1970),  we  assume  that  a  power  function  relates  ground  pressure  to  sinkage.  The 
SSMl  .0  calculates  sinkage  (referenced  to  the  original  snow  depth  hj  for  a  given  station  from 

2  =2final  (AS) 

(see  Mellor  1964,  Fig.  III-34). 

The  sinkage  z.,  at  a  given  station  /  on  the  vehicle,  can  then  be  calculated  from  (Fig.  A3) 


2i  -  Zfinal  (pi/pmax)  ~  Z 

i-1 

Zi  =  0 


((Pi  >  Pi-1.  Pi-2  •  ■  •) 
(Pi  ^  Pi-l.  Pi-2---) 


We  now  have  all  of  the  parameters  necessary  to  calculate  T  using  eq  A3  (i.e.,  W,  p  [or  f],  b, 
r,  z,  c  and  (})). 

Recognizing  that  the  snow  properties  are  different  in  the  flat  and  curved  portions  of  the 
contact  zone,  we  divide  the  traction  equation  into  two  parts 

Tj  =AijCi_]  +A2:Cj  + — !- Wj  tan<t)i_i  +  — ZL  W;  tanc));  (AlO) 

Aj  Ai 

where  c.  ^  and  (fij  j  are  the  snow  strength  parameters  for  the  snow  in  the  curved  contact  region 
A 1  ,  and  c.  and  ()).  are  for  the  snow  directly  under  the  tire  or  track  at  station  ;  (Fig.  A3). 
Equations  A4  are  used  with  p.  j  to  obtain  c.  j  and  p.  j,  and  with  p.  to  get  c.  and  (])..  This 
expression  assumes  that  the  normal  load  supported  by  A  i  is  proportional  to  the  ratio  of 
that  area  to  the  total  area  A..  The  load  on  A 2  is  treated  similarly. 

To  calculate  motion  resistance  we  need  to  know  the  intermediate  values  of  snow 
density  as  compaction  progresses  from  initial  density  p^  to  pj.  We  have  already  stated  that 
Pjis  associated  with  maximum  sinkage  z  =  z^^^and  p^  corresponds  to  a  sinkage  of  z  ==  0. 
Recalling  eq  A7,  and  the  assumptions  that  it  is  based  on,  we  can  find  the  density  beneath  a 
particular  station  /  from 
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(All) 


Pi  = - - 

1  -(zj//;  i_i) 

where  pj  j  and  It._^  are  the  snow  density  and  depth  prior  to  the  passage  of  the  tire  or  track  at 
station  i,  and  z.  (eq  A9)  is  the  sinkage  produced  at  this  station  (Fig.  A3). 

Lastly,  we  recognize  that  not  all  of  the  tires  or  tracks  on  a  particular  vehicle  may  be 
traveling  in  undisturbed  snow.  Some  tires  or  tracks  may  follow  exactly  in  the  path  of  a 
preceding  element,  or  may  operate  in  undisturbed  snow  (e.g.,  a  narrow  or  wide  trailer 
behind  a  vehicle)  or  may  encounter  both  compacted  and  uncompacted  snow  (e.g.,  dual  tires 
following  a  single  tire).  We  need  then  to  account  for  the  possibility  of  tires  (tracks)  having 
some  percentage  (a)  of  their  width  compacting  new  snow  while  the  remainder  is  traveling 
in  a  previously  created  rut.  The  equations  for  A,,  R  and  T  are  therefore  modified  to  become 

IS  g 


H - 1 - H 
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Aij  =a 


bi  rj  arccos 


(  i  \ 

n-1 2j 


\  n  /j 


(l-o) 


bi  r,  arccos  I  ^ - - 

f  i 


(A12) 


Rsi=a 


Pibiho  Po(— J —  In^  -i’ 
Pi  -po  Po  Pi 


+  (1-a) 


pbi  h  i_i  pi_i  — 1 - In^  -  i 

Ipi-Pi-I  Pi-1  Pi 


Pi  >  Pi-1 


Rs-  =a 

PihhoPoi— — 

In  Pi- 

Ipi-Po 

Po 

Pi 

Pi  —  Pi-l  ^  Po 


Rsi  =  0 


Pi  —  Po 


(A13) 


A-y- 

Tgj  =A2jCi  +1 — !-Wi|tan(l)i  +a 

\  Ai 


^\iCo  +(’— Wi|  tan())o 
^i 


+  (1-a; 


^li  ^i-1  Wi|tan(t)i_i 

Ai 


(A14) 


for  each  station.  Equations  A13  and  A 14  provide  the  essence  of  the  SSMl.O.  These  equations  are 
executed  foreach  station  of  the  vehicle  and  a  running  sum  for  traction  and  resistance  accumulated.  The 
ultimate  ability  of  the  vehicle  to  move  (total  net  traction,  is  then  determined  from 

Tn  =  l(Tg.-R^)  (A15) 

i=l 

where  n  is  the  number  of  tires  or  tracks  on  the  vehicle.  If  T^is  positive,  the  vehicle  is  mobile 
and  has  the  capacity  to  accelerate,  climb  hills  or  pull  a  payload  in  proportion  with  the 
magnitude  of  T^.  A  value  of  zero  indicates  impending  immobilization,  and  a  negative  value 
of  predicts  a  definite  no-go  situation.  A  copy  of  the  SSMl.O,  in  HP  Basic  computer  code 
is  contained  in  the  next  section,  along  with  the  output  from  sample  runs,  and  a  flow  chart  is 
provided  in  Figure  A4. 
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COMPUTER  CODE  AND  SAMPLE  OUTPUT 


10  I 

20  I  SSM/1 .01 

•  (trig  tr  slnkage  and  snow  density  function  of  cont .  pr.) 

40  !  (manual  input  of  all  vehicle  and  snow  data) 

B0  I  (Liston’s  resis;  best-fit  c,  phi,  tuo-area  traction) 

£0  .1 

70  net“0 
80  RR-0 
90  DB-0 
100  go'0 

110  p 1 «0 

120  z»0 
130  sumz=0 

140  DISP  "vehicle'^" 

1S0  INPUT  UehS 

160  DISP  "wheeled  (w)  or  tracked  (t)"?" 

170  INPUT  types 

180  DISP  "highest  station  ground  pressure  (kPa)«" 

190  INPUT  pmax 

200  pmax«pmax»0 . 1  i  convert  from  kPa  to  N/cm“2 
210  I 

220  I  ••  establish  final  density  based  on  largest  footprint  pressure  •» 

230  sigmaf“0.5 

240  IF  pmax>21  THEN  5igmaf»0.55 
250  IF  pmax>35  THEN  sigmaf“0.6 
260  IF  pmax>70  THEN  8igmaf«0.65 
270  ! 

280  DISP  "snow  depth  (cm)  «“ 

290  INPUT  h0 

300  DISP  "snow  density  (g/cm"3)  «* 

310  INPUT  signa0 
320  sigmal  ••sigma0 
330  signa2"sigms0 
340  h1mh0 
3S0  I 

360  PRINT  "Vehicle!  "iVehS 
370  PRINT  *  ‘ 

380  PRINT  USING  formatl  i  "initial  state:  snow  depth  ""ih0i"  cm” 

390  PRINT  USING  format2  i  "  snow  density  •’i5igma0i"  g/cm*3" 

400  ! 

410  zmax>h0*( 1 -sigma0/sigmaf  )  I  calculate  maximum  sinkage 
420  IF  zmax<0  THEN  zmax«0 
430  ! 

440  1  ••  calculate  snow’s  cohesion  based  on  the  Initial  snow  density  •» 

450  coh0“0 , 666+88 . 7»s igma0 

460  coh0”coh0» 1 000  I  convert  from  kPa  to  Pa 

470  IF  coh0>27270  THEN  coh0«27270  I  clip  high  cohesion  values 

480  I 

490  1  ••  calculate  snow’s  angle  of  friction  based  on  initial  density  •• 

500  phi0«17. 06-16. 3»signa0 

510  phi0«phi0»(PI/180 )  I  convert  from  degrees  to  radians 

520  ! 

530  coh1»coh0 
540  phi1“phi0 
550  coh2”coh0 
560  phi2~phi0 
570  I 

580  RAD  !  compute  in  radians  units  for  trig  functions 
590  I 

600  PRINT  ■  " 

610  PRINT  USING  format!  t  "  cohesion  -  "icoh0("  N/m*2  (Pa)" 

620  PRINT  USING  format2  i  "  angle  of  friction  -  "iphi0i"  radians 

630  PRINT  "  ■ 

640  I 

650  !  ••  enter  vehicle  data  one  tire  station  at  a  time  •• 
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B60  OISP  "how  many  wheel  or  track  stations  on  each  side  of  the  vehicle?" 

670  input  N 
680  I 

690  FOR  I«1  TO  N 
700  OISP  ““ 

710  OISP  "station  " t I 

720  OISP  “wheel  radius  or  approx  radius  of  compacting  portion  of  track  cm)  »“ 
730  INPUT  rads 

740  OISP  "single  tire  or  track  width  at  this  location  (cm) 

750  INPUT  wid 

760  OISP  "single  tire  or  track  load  <N)  •' 

770  INPUT  loa 

780  IF  type$-"t"  THEN  GOTO  skipl 
790  OISP  “inflation  pressure  ( kPa  ) 

800  INPUT  pres 

810  pres“pres*0 . 1  '  convert  from  kPa  to  N/cm''2 
820  GOTO  skip2 
830  I 

840  skipl:  OISP  "length  of  flat  portion  of  track  (cm)" 

850  INPUT  tl 

860  pres“loa/(  1 1 ‘wid  )  i  calculate  ground  pressure  in  N/cm''2 
870  I 

880  skip2:  OISP  "driven?" 

890  INPUT  rS 

900  IF  type$-"t"  THEN  GOTO  skipS 
910  OISP  "duals  ?• 

920  INPUT  Rt 

930  IF  R$<>"y"  AND  RS<>"Y"  THEN  SOTO  970 
940  utd“wid«2 
950  loa“loa»2 
960  I 

970  8kip3:  OISP  "percent  of  width  compacting  virgin  snow  <X) 

980  INPUT  prcnt 
990  I 

1000  IF  pres>pl  AND  sigma2<sigmaf  THEN  GOTO  1060  !  check  for  added  sink,  here 

1010  z-0 

1020  tenpres“0 

1030  IF  prcnt>0  AND  sigma2>8igma0  THEN  GOTO  1180 
1040  re5is"0 
1050  GOTO  1240 

1060  z“(pres/pmax  )“0.5»znax-sumz  1  calculate  additional  sinkage  this  station 
1070  I 

1080  I  set  rut  bottom  values  •• 

1090  h2”hl-z 

1100  slgma2”signal /(  1-z/hl  ) 

1110  phi2-( 17.06-16.3»sigma2)*(PI/l80) 

1120  coh2”<0. 666+88 . 7» s igma2  )• 1 000 

1130  IF  coh2>27270  THEN  coh2“27270  >  clip  high  cohesion  values 
1140  I 

1150  I  ••  calculate  resistance  parameter  •• 

1160  tempres“l /< 5 igma2-sigmal  )»L0G( sigma2/sigma1  )- 1 /s igma2  '  in  rut 
1170  I 

1180  I  ••  calculate  motion  resistance  at  this  station  •• 

1190  res i s*prcnt / 1 00» ( pres»wi d*h0* 5 igma0 )  •  in  virgin  snow 

1200  resls»resis  +  < 1 /< sigma2-sigma0  )»L0G( signa2/sigma0 )-l /5igma2  )  'virgin  snow 
1210  moreres" tempres* ( 1 -prcnt / 1 00 )• ( pres»wi d»h 1 *51 gma I  )  '  in  rut 
1220  resis-resis+moreres 
1230  I 

1240  I  calculate  area  in  contact  with  the  snow  ♦♦ 

1250  area  1 "prcnt / 1 00* ( wi d»rads» ACS( ( rads-(  sumz  +  z ) )/rads  )  )  '  In  virgin  snow 

1260  area  1 “area  1 +( 1 -prcnt / 1 00  )•< wi d»rads»ACS( ( rads-z ) /rads  )  )  '  in  cm'2 

1270  area  1 “area I / 1 00" 2  '  convert  to  m"2 

1280  aroa2“loa/pres  !  in  cm"2 

1290  eroa2“area2/ 1 00"2  I  convert  to  m"2 

1300  area“areal +area2 

1310  I 

1320  '  ••  calculate  gross  traction  at  this  station  •• 

1330  trac“area I ‘cohl +area2»coh2 
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1 340  trac-trac+area 1 /area* loa»TAN( phi  1  )  +  arca2/area* loa*TAN( phi2  > 

1350  IF  type$"“t"  THEN  GOTO  1380  •  negate  firat  time  adjustment 
1 3G0  IF  I"I  THEN  trac“area*coh1 +Joa*TAN(phi 1  )  !  no  rut  bottom  adjustment 
1370  ! 

1380  iF  r$“"n"  OR  r$«"N"  THEN  trac*0  !  no  traction  if  not  driven 
1390  I 

1400  !  ••  double  to  account  for  both  sides  of  the  vehicle  *• 

1410  area“2*area 
1420  trac»2»trac 
1430  re8i5“2*re8i8 

1440  ! 

1450  1  ••  print  output  for  this  station  •• 

1460  a1“2*area1  I  for  printout 
1470  a2”2*area2  !  for  printout 
1480  PRINT  "  " 

1490  PRINT  "station  "il 

1500  PRINT  USING  formatl  i  “  additional  sinkage  ■"izi"  cm" 

1510  PRINT  USING  format2  \  '  total  area"" i area » "  m*2  (curved  ""lalt"  .flat  “ " 
ia2i"  )• 

1520  PRINT  USING  formatl  i  "  snoui  resis  “"iresisi"  N" 

1530  PRINT  USING  formatl  t  ’  snou  trac  -"itraci"  N" 

1540  I 

1550  I  ••  running  summation  for  whole  vehicle  •• 

1560  RR"RR+resi5  '  sum  for  whole  vehicle 
1570  DB"DB+trac  I  sum  for  whole  vehicle 

1580  net“trac-resis  i  calculate  net  traction  for  individual  station 
1590  PRINT  USING  formatl  !  "  net  snow  traction  for  station  *  "met;"  N" 

1600  PRINT  USING  formatl  i  "  rut  bottom:  depth*" ih2;"  cm" 

1610  PRINT  USING  format2  i  "  densi t y " t s igna2 t ’  g/cm'3" 

1620  PRINT  USING  formatl  i  "  cohes ion» " i coh2 ; ”  N/n“2" 

1630  PRINT  USING  format2  i  "  f r i c t ion* " i phi 2  i "  radians" 

1640  Qo*go+r,et  I  sum  net  traction  for  vehicle 

1650  > 

1660  I  ♦»  save  last  stat  on  values  for  next  iteration  •• 

1670  IF  pres>pl  THEN  pl*pres 

1680  sigmal *5igma2 

1690  hl-h2 

1700  sumz=sum2+z 

1710  cohl=coh2 

1720  phi1«phi2 

1730  I 

1740  NEXT  I 

1750  I 

1760  I  ••  calculate  mobility  in  English  units  for  output  •• 

1770  eRR-RR/4. 448222 
1780  eDB"DB/4. 448222 
1790  ego*go/4. 448222 
1800  ez*sumz/2.54 
1810  I 

1820  I  print  out  whole  vehicle  results  •• 

1830  PRINT  ■  " 

1840  PRINT  "  • 

1850  PRINT  USING  formatl  i  "  total  sinkage  for  vehicle*  "isumzi"  cm  ("lezi"  in 

)” 

1860  PRINT  USING  formatl  t  "  total  snow  resistance  *  “tRRt"  N  ("leRR;"  lb)" 
1870  PRINT  USING  formatl  i  "  total  snou  traction  “’iDBi"  N  (“leDBi"  lb>" 

1880  PRINT  '  ■ 

1890  PRINT  USING  formatl  i  "  net  traction  for  vehicle  *  “igot"  N  ("tegoi"  lb)” 
1900  PRINT  ■  " 

1910  PRINT  •  • 

1920  PRINT  •  ■ 

197'*  ‘ 

1940  formatl:  IMAGE  3< K .DDODOO.OO  ) 

1950  format2:  IMAGE  3<K.0Z.DDDD) 

I960  I 
1970  END 
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Uehicla:  HMMUU 


initial  state:  snou  depth  “  5.00  cm 

snow  density  «  0.5500  Q/cn“3 

cohesion  ■  27270.00  N/n^2  (Pa) 

angle  of  friction  •  0.1413  radians 


station  1 

additional  sinkage  ”  0.00  cn 

total  area”  0.1048  n“2  (curved  “  0.0000  .flat  “  0 
snou  resis  •  0.00  N 

snou  trac  ”  4915.01  N 

net  snou  traction  for  station  •  4915.01  N 

rut  bottom:  depth-  0.00  cm 

density-  0.5500  g/cm''3 
cohesion*  27270.00  N/m''2 
friction-  0.1413  radians 


station  2 

additional  sinkage  -  0.00  cm 

total  area-  0.1252  m“2  (curved  *  0.0000  .flat  *  0 
snou  resis  *  0.00  N 

snou  trac  *  6115.90  N 

net  snou  traction  for  station  *  6115.90  N 

rut  bottom:  depth-  0.00  cm 

density-  0.5500  g/cm"3 
cohesion-  27270.00  N/m''2 
friction-  0.1413  radians 


total  sinkage  for  vehicle-  0.00  cm  (  0.00 
total  snou  resistance  •  0.00  N  (  0.00  lb) 
total  snou  traction  •  11030,92  N  (  2479.85  lb) 

net  traction  for  vehicle  -  11030.92  N  (  2479.85 


.  1048  ) 


.  1252  ) 


in  ) 


lb  ) 
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Vehicle:  HEMTT  <241/276  > 


initial  state:  snow  depth  -  25.00  cm 

snow  density  •  0.2200  Q/cn'3 

cohesion  »  20180.00  N/i<i“2  (Pa) 

angle  of  friction  »  0.2352  radians 

station  1 

additional  sinkage  -  14.03  cn 

total  area=  0.6517  n“2  (curved  =  0.3947  .flat  =  0.2570) 
snow  resis  «  11515.55  N 
snow  trac  »  28007.95  N 
net  snow  traction  for  station  “  16492.41  N 

rut  bottom:  depth*  10.97  cm 

density*  0.5014  g/cm"3 
cohesion*  27270.00  N/n'2 
friction*  0.1551  radians 

station  2 

additional  sinkage  *  0.00  cm 

total  area*  0.2608  m*2  (curved  -  0.0000  .flat  -  0.2608) 
snow  resis  •  0.00  N 

snow  trac  «  16949.49  N 
net  snow  traction  for  station  *  16949.49  N 

rut  bottom:  depth*  10.97  cm 

density*  0.5014  g/cm*3 
cohesion*  27270.00  N/m“2 
friction*  0.I55I  radians 

station  3 

additional  sinkage  *  .97  cm 

total  area-  0.3636  m'2  (curved  -  0.1040  .flat  -  0.2595) 
snow  resis  *  1230.94  N 

snow  trac  -  20383.43  N 
not  snow  traction  for  station  -  19152.48  N 

rut  bottom:  depth*  10.00  cm 

density*  0.5500  g/cm“3 
cohesion*  27270.00  N/m''2 
friction*  0.1413  radians 

station  4 

additional  sinkage  *  0.00  cm 

total  area*  0.2613  m“2  (curved  *  0.0000  .flat  *  0.2613) 
snow  resis  *  0.00  N 

enow  trac  *  17374.44  N 
net  snow  traction  for  station  «  17374.44  N 

rut  bottom:  depth-  10.00  cm 

density*  0.5500  g/cm“3 
cohesion*  27270.00  N/m"2 
friction*  0.1413  radians 


total  sinkage  for  vehicle*  15.00  cm  (  5.91  in) 

total  snow  resistance  *  12746.49  N  (  2865.52  lb) 

total  snow  traction  *  82715.30  N  (  18595.14  lb) 

net  traction  for  vehicle  *  69968.82  N  (  15729.61  lb) 
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Vehicle:  SUSV 


initial  state:  snow  depth  «  9.00  cn 

snou  density  •  0.1480  a/cn"3 

cohesion  »  13793.60  N/n"2  (Pa) 

angle  of  friction  -  0.25B6  radians 


station  1 

additional  sinkage  “  6.35  cn 

total  area=  2.5505  n"2  (curved  «  0.2280  .flat  *  2.3226) 
snoui  resis  •  313.84  N 

snou  trac  =  71574.85  N 
net  snou  traction  for  station  “  71261.01  N 

rut  bottom:  depth”  2.65  cn 

density”  0.5024  g/cn''3 
cohesion”  27270.00  N/n''2 
friction”  0.1548  radians 


stat  ion  2 

additional  sinkage  ■  0.00  cn 

total  area-  2.3226  n''2  (curved  •  0.0000  .flat  -  2.3226) 
snou  reals  •  0.00  N 

snow  trac  *  68106.12  N 
net  snou  traction  for  station  =  68106.12  N 

rut  bottom:  depth-  2.65  cm 

density”  0.5024  Q/cn*3 
cohesion-  27270.00  N/n‘2 
friction-  0.1548  radians 


total  sinkage  for  vehicle-  6.35  cm  (  2.50  in) 

total  snou  resistance  •  313.84  N  (  70.55  lb) 

total  snou  traction  -139680.97  N  (  31401.53  lb) 

net  traction  for  vehicle  -  139367.14  N  (  31330.98  lb) 
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APPENDIX  B:  TEST  VEHICLE  DATA 


CIV,  179  kPa  WHEELED 
GVW  -  29696.4  H 

MAXIMUM  GROUND  (INFLATION)  PRESSURE  179.30  kP« 


AVERAGE 

HARD  SURFACE 

CONTACT  AREA 

0.026 

WEIGHT  PER 

TIRE 

NEW  SNOW 

STATION 

RADIUS  (co) 

WIDTH  (cm) 

WHEEL  (N) 

PRES.  (kPa) 

DRIVEN 

duals 

PERCENT 

1 

35.6 

25.4 

6174.1 

179.3 

Y 

N 

100 

2 

35.6 

25.4 

6174.1 

179.3 

N 

N 

0 

CIV,  103  kPa  WHEELED 

GVW  -  : 

24696.4  N 

MAXIHUM 

GROUND  (INFLATION)  PRESSURE  110.30  kPa 

AVERAGE 

HARD  SURFACE 

CONTACT  AREA 

0.0345  B* 

WEIGHT  PER 

TIRE 

NEW  SNOW 

STATION 

RADIUS  (cb) 

WIDTH  (cm) 

WHEEL  (N) 

PRES.  (kPa) 

DRIVEN 

DUALS 

PERCENT 

1 

34.9 

26.7 

6174.1 

110.3 

Y 

N 

100 

2 

34.9 

26.7 

6174.1 

110.3 

N 

N 

0 

HEMTT,  207/276  kPa  WHEELED 

GW  -  268560.1  N 

MAXIMUM 

GROUND  (INFLATION)  PRESSURE  275.80  kP« 

AVERAGE 

HARD  SURFACE 

CONTACT  AREA 

0.149 

WEIGHT  PER 

TIRE 

NEW  SNOW 

STATION 

RADIUS  (cb) 

WIDTH  (CB) 

WHEEL  (N) 

PRES.  (kPa) 

DRIVEN 

DUALS 

PERCENT 

1 

61.7 

46.5 

31004.0 

241.3 

Y 

N 

100 

2 

61.7 

46.5 

31459.9 

241.3 

Y 

N 

0 

3 

61.7 

47.5 

35785.8 

275.8 

Y 

N 

0 

4 

61.7 

47.5 

36030.4 

275.8 

Y 

N 

0 

HEMTT.  138/207  kPa  WHEELED 

GVW  -  268560.1  N 

MAXIMUM 

GROUND  (INFLATION)  PRESSURE  206.80  kPa 

AVERAGE 

HARD  SURFACE 

CONTACT  AREA 

0.171 

WEIGHT  PER 

TIRE 

NEW  SNOW 

STATION 

RADIUS  (cb) 

WIDTH  (cb) 

WHEEL  (N) 

PRES.  (kPa) 

DRIVEN 

DUALS 

PERCENT 

1 

58.9 

48.9 

31004.0 

137.9 

Y 

N 

100 

2 

58.9 

48.9 

31459.9 

137.9 

Y 

N 

0 

3 

58.9 

48.3 

35785.8 

206.8 

Y 

N 

0 

4 

58.9 

48.3 

36030.4 

206.8 

Y 

N 

0 

HMMWV,  138/152  kPa  WHEELED 

GW  -  ; 

13450.5  N 

MAXIMUM 

GROUND  (INFLATION)  PRESSURE  151.70  kFa 

AVERAGE 

HARD  SURFACE 

CONTACT  AREA 

0.074  ■» 

WEIGHT  PER 

TIRE 

MEW  SNOW 

STATION 

RADIUS  (ca) 

WIDTH  (cb) 

WHEEL  <N) 

PRES.  (kPa) 

DRIVEN 

DUALS 

PERCENT 

1 

42.9 

32.0 

7228.3 

137.9 

Y 

N 

100 

2 

42.9 

33.0 

9496.9 

151.7 

Y 

N 

0 

LAV.(12.SX20)  207  kPa  WHEELED 
GVW  -  125483.7  N 

MAXIMUM  GROUND  (INFUTION)  PRESSURE  206.80  kPa 
AVERAGE  HARO  SURFACE  CONTACT  AREA  0.100 
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STATION 

RADIUS  (ca) 

WIDTH  (ca) 

WEIGHT  PER 
WHEEL  (N) 

TIRE 

PRES.  (kPa) 

DRIVEN 

DUALS 

NEW  SNOW 
PERCENT 

1 

44.5 

35.3 

16102.5 

206.8 

Y 

N 

100 

2 

44.5 

35.3 

17281.3 

206.8 

Y 

N 

0 

3 

44.5 

35.3 

14501 . 1 

206.8 

Y 

N 

0 

4 

44.5 

35.3 

14857.0 

206.8 

y 

N 

0 

LAV.  <12 

5X20)  103  kPa  WHEELED 

GVW  -  125483.7  N 

MAXIMUM 

GROUND  (INFLATION)  PRESSURE  103.40  kPa 

AVERAGE 

HARD  SURFACE 

CONTACT  AREA 

0.141  a> 

WEIGHT  PER 

TIRE 

NEV  SNOW 

STATION 

RADIUS  (ca) 

WIDTH  (ca) 

WHEEL  (H) 

PRES.  (kPa) 

DRIVEN 

DUALS 

PERCENT 

1 

41.4 

37.8 

16102.5 

103.4 

Y 

N 

100 

2 

41.4 

37.8 

17281.3 

103.4 

Y 

N 

0 

3 

41.4 

37.8 

14501.1 

103.4 

Y 

N 

0 

4 

41.4 

37.8 

14857.0 

103.4 

Y 

N 

0 

LAV, (11X16)  289  kPa 

WHEELED 

GVU  -  119634.3  N 

MAXIMUM 

GROUND  (INFLATION)  PRESSURE  289.60  kPa 

AVERAGE 

HARD  SURFACE 

CONTACT  AREA 

0.058  a> 

HEIGHT  PER 

TIRE 

NEW  SNOW 

STATION 

RADIUS  (ca) 

WIDTH  (ca) 

WHEEL  (N) 

PRES.  (kPa) 

DRIVEN 

DUALS 

PERCENT 

1 

43.4 

31.4 

15346.3 

289.6 

Y 

N 

100 

2 

43.4 

31.4 

16480.6 

289.6 

Y 

N 

0 

3 

43.4 

31.4 

13822.8 

289  6 

Y 

N 

0 

4 

43.4 

31.4 

14167.5 

289.6 

Y 

N 

0 

LAV, (11X16)  165  kPa 

WHEELED 

CVH  -  119634.3  N 

MAXIMUM 

GROUND  (INFUTION)  PRESSURE  165.50  kPa 

AVERAGE 

HARD  SURFACE 

CONTACT  AREA 

0.102  8> 

WEIGHT  PER 

TIRE 

NEW  SNOW 

STATION 

RADIUS  (ca) 

WIDTH  (cm) 

WHEEL  (N> 

PRES.  (kPa) 

DRIVEN 

DUALS 

PERCENT 

1 

41.7 

33.2 

15346.3 

165.5 

Y 

N 

100 

2 

41.7 

33.2 

16480.6 

165.5 

Y 

N 

0 

3 

41.7 

33.2 

13822.8 

165.5 

Y 

N 

0 

4 

41.7 

33.2 

14167.5 

165.5 

Y 

N 

0 

SUSV  TRACKED 
GVW  -  61340.7  N 

MAXIMUM  GROUND  (INFLATION)  PRESSURE  13.20  kPa 
AVERAGE  HARD  SURFACE  CONTACT  AREA  1.180  ■’ 


VfEIGHT  PER  TRACK  NEV  SNOW 


STATION 

RADIUS  (ca) 

WIDTH  (ca) 

TRACK  (N)  LENGTH 

(ca) 

PERCENT 

1 

26.40 

60.96 

15390.80 

190.50 

100.00 

2 

26.40 

60.96 

15279.60 

190.50 

0.00 

H113  TRACKED 

CVH  -  104087.9  N 

MAXIMUM  GROUND  (INFLATION) 

PRESSURE  51.02  kPa 

AVERAGE  HARD  SURFACE  CONTACT  AREA  1.020  ■* 

HEIGHT  PER  TRACK  NEW  SNOW 

STATION  RADIUS  (ca)  HIDTH  (ca)  TRACK  (N)  LENGTH  (ca)  PERCENT 

1  36.80  38.10  32043. SO  266.70  100.00 
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BRADLEY  TRACKED 
GW  -  223299.6  N 

MAXIMUM  GROUND  (INFLATION)  PRESSURE  53.09  kP* 

AVERAGE  HARD  SURFACE  CONTACT  AREA  2.090 

HEIGHT  PER  TRACK  NEV  SNOW 

STATION  RADIUS  (ca)  HIDTH  (cm)  TRACK  (N)  LENGTH  (ca)  PERCENT 


1  35.60  53.39  111649.75  391.20  100.00 


5 -TON,  207  kP«  WHEELED 
GW  -  105778.2  N 

MAXIMUM  GROUND  (INFLATION)  PRESSURE  206.80  kPa 
AVERAGE  HARD  SURFACE  CONTACT  AREA  0.171  a> 


HEIGHT  PER  TIRE  NEW  SNOW 

STATION  RADIUS  (ca)  HIDTH  (ca)  WHEEL  (N)  PRES.  (kPa)  DRIVEN  DUALS  PERCENT 


1  55.4  40.8 

2  55.4  39.2 

3  55.4  39.2 


23842.4 

14523.4 
14523.4 


206.8 

206.8 

206.8 


Y  N  100 

Y  N  0 

Y  N  0 
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Drawbar  F’ull  Coefficient  S  Drawbar  Pull  Coefficient  Drawbar  Pull  Coefficient 
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APPENDIX  D:  SNOW  DATA 


UNDISTURBED  SNOW  CHARACTERIZATION 


DATE _ Qi-21-88 _ 

TIME  _ 

AMB  TEMP  (^'O  _ 

SITE  NAME  SLOPES 


DEPTH  (Cl) 

CRYSTAL  SIZE  (mo) 

ICE  THICKNESS  lea) 

10-56 

1-1.5 

DEPTH  (cal 

DENSITY  (9/c«^ 

TEHPEHATURE  (®C1 

10 

0.09 

-3 

20 

0.26 

-1 

30 

0.30 

0 

40 

0.28 

0 

56 

ground 

Note:  Depths  are  measured  down  from  the  pack  surface 

DATE _ lrlZ-8a _ 

TTMF  15:00 _ 

AMB  TEMP  (°C)  -5-5 _ 

SITE  NAME  tlnk _ 


DEPTH  (Cl) 

CRYSTAL  SIZE  (an) 

ICE  THIOOCSS  (ca) 

0-14 

1  fllahcT 

DEPTH  (cal 

DENSITY  (e/ei^ 

TEHPEHATURE  (“0 

0-3 

<0.06 

5 

-5 

10-13 

0.10 

14 

arounil 

Note:  Depths  are  measured  down  from  the  pack  surface 
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UNDISTURBED  SNOW  CHARACTERIZATION 

nATF  01-23-88 _ 

TIME _ lii” _ 

AMB  TEMP  (°C)  --a.  2 _ 

SITE  NAME  _ 


Note:  Depths  are  measured  down  from  the  pack  surface 


DATE  i-26-gB _ 

TTMF  1^:00 _ 

AMB  TEMP  (Oc)  -13 
SITE  NAME  Runway _ 


DEPTH  (Cl)  DTCITY  (g/ca^)  TEHPiSATUHE  (®CI  DEPTH  Ici)  CBYSTM.  SIZE  W  ICE  THICKNESS  (ca) 


Note;  Depths  are  measured  down  from  the  pack  surface 
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UNDISTURBED  SNOW  CHARACTERIZATION 


Note;  Depths  are  measured  down  from  the  pack  surface 

HATF  02-10-88 _ 

TIME  12  =  0° _ 

AMB  TEMP  (OC)  -17-0 _ 

SITE  NAME  1 _ 


Note:  Depths  are  measured  down  from  the  pack  surface 


UNDISTURBED  SNOW  CHARACTERIZATION 


QAJE  02-10-88 _ 

TTMF  1^  =  00 _ 

AMB  TEMP  (°C)  -lA _ 

SITE  NAME  rink  2  (Entrance  End) 


Note:  Depths  are  measured  down  from  the  pack  surface 


AMB  TEMP  (OC)  -1* 
SITE  NAME  RINK  2  Cent«r 


Note:  Depths  are  measured  down  from  the  pack  surface 
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UNDISTURBED  SNOW  CHARACTERIZATION 


DATE _ _ 

TIME _ lilQfl _ 

AMB  TEMP  (°C)  _-l» 
SITE  NAME  rink  2 


DEPTH  (Cl)  CRYSTAL  SIZE  Ilia)  ICE  THICKNESS  lea) 


Note:  Depths  are  measured  down  from  the  pac)t  surface 

DATE  02-11-88 

TTMF  16:15 _ 

AMB  TEMP  (°C)  jri* _ 

SITE  NAME  ^ _ 


DEPTH  In)  CRYSTAL  SIZE  laa)  ICE  THICKIESS  lea) 


Note:  Depths  are  measured  down  from  the  pack  surface 
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UNDISTURBED  SNOW  CHARACTERIZATION 


DATE _ 2-12-88 _ 

TTMF  13:10 _ 

AMB  TEMP  (°C)  -16 

SITE  NAME  _ Rink  2 


DEPTH  (c«) 

DENSITT  tg/c«3) 

1-1, 

0.12 

Note:  Depths  are  measured  down  from  the  pack  surface 

□  ATF  02-12-88 _ 

TIMF  13:20 _ 

AMB  TEMP  (°C)  -16 _ 

SITE  NAME  _ Rink  1 _ 


Note:  Depths  are  measured  down  from  the  pack  surface 
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UNDISTURBED  SNOW  CHARACTERIZATION 


QAJE  02-13-88 _ 

TTMF  9:50 _ 

AMB  TEMP  (°C) 

SITE  NAME  RIULl— 


nATF  02-13-88 _ 

TTMF  10:5t _ 

AMB  TEMP  (°C) 
SITE  NAME  RINK  2 


Note:  Depths  are  measured  down  from  the  pack  surface 
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UNDISTURBED  SNOW  CHARACTERIZATION 


DATE _ o^-U-g8 

TTMF  11:00 _ 

AMB  TEMP  (°C)  ^ 
SITE  NAME  SLOPES 


DEPTH  (ce) 

DENSITY  |g/c«3) 

TEMPEHATUBE  (®CI 

<5-1? 

0.22 

22-25 

0.35 

40-43 

0.36 

57-60 

0.34 

83 

ground 

DEPTH  tee) 

CRYSTAL  SI2E  (en) 

ICE  THICKNESS  (ce) 

Note:  Depths  are  measured  down  from  the  pack  surface 

nATF  02-13-86 _ 

TTMF  llria _ 

AMB  TEMP  (°C)  _ 

SITE  NAME  -  Ttxag  cle»r«d  «rt« _ 


DEPTH  tea) 

CRYSTAL  SIZE  M 

ICE  THICKNESS  tee) 

0-23 

1  ftiehc) 

23 

1 

24-39 

1-1.5  (looie) 

DEPTH  tea) 

DENSITY  (g/ea^ 

TEW>£HATUHE  (*« 

1-3 

0. 14 

8-11 

0.26 

10 

-13 

17-20 

0.27 

25 

-6 

26-29 

0.23 

31-34 

0.25 

32 

-2 

39 

ground 

_ 1 

Note:  Depths  are  measured  down  from  the  pack  surface 
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UNDISTURBED  SNOW  CHARACTERIZATION 


DATE _ QZ-n-88 _ 

TTMP  11:18 _ 

AMB  TEMP  (°C)  ^ _ 

SITE  NAME  Texas  uncleared 


DEPTH  (c«)  DENSITY  (g/cB^  TEMPEHATUPE  (“Cl  j  DEPTH  (ce)  CBYSTAL  SIZE  (lu)  ICE  THICKNESS  (ca) 


Note;  Depths  are  measured  down  from  the  pack  surface 


nATP  02-13-88 _ 

TIME _ JlHi - 

AMB  TEMP  (OC)  - 

SITE  NAME  South  End  Runvay 


Note:  Depths  are  measured  down  from  the  pack  surface 
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UNDISTURBED  SNOW 


raTF  02-13-88 _ 

TIME _ iliiZ _ 

AMB  TEMP  (°C)  -10 

SITE  NAME  North  End  Runway 


DEPTH  (ca) 

DENSITY  l9/c»^ 

TENPEHATURE  (“Cl 

0-3 

0.20 

5-8 

0.29 

5 

-13 

13-16 

0.20 

14 

-10 

17-20 

0.26 

20 

-8 

26 

eround 

Note:  Depths  ere  measured  down  from  the  pack  surface 


DATE _ 02-15-88 _ 

TTMF  K:** _ 

AMB  TEMP  (OC)  -11.0 
SITE  NAME  MNK  2 


DEPTH  lea) 

DENSITY  tg/ca^ 

TEI«>EMTURE  (“O 

0-3 

0.20 

4-7 

0.26 

5 

-7 

8-11 

0.29 

12-15 

0.27 

15 

-5 

19-22 

0.23 

26-27 

0.26 

25 

-3 

28-31 

0.28 

31 

ground 

Note:  Depths  are  measured  down  from  the  pack  surface 
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ZATION 


UNDISTURBED  SNOW  CHARACTERIZATION 


riATF  Q2-iR-aR _ 

TIME _ LiifiS _ 

AMB  TEMP  (Oc)  *6.5 
SITE  NAME  RIUK  1 


Note:  Depths  are  measured  down  from  the  pack  surface 
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UNDISTURBED  SNOW  CHARACTERIZATION 


date  02-16-88 _ 

TTMF  10  =  00 _ 

AMB  TEMP  (°C)  -  -1— 

SITE  NAME  N««r  Slopea 


Note;  Depths  are  measured  down  from  the  pack  surface 


Note:  Depths  are  measured  down  from  the  pack  surface 
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UNDISTURBED  SNOW  CHARACTERIZATION 


DATE _ 02-18-1988 

TTMF  11:90 _ 

AMB  TEMP  (°C)  _±I 
SITE  NAME  mwK  7 


DEPTH  (cb) 

CRYSTAL  SIZE  (aal 

ICE  THICKNESS  lea) 

0-16 

1  (tight) 

16 

1 

17-28 

1-2  (loose) 

DEPTH  lea) 

DENSITY  (5/cb3| 

TEMPERATURE  (“Cl 

0-3 

0.26 

4-7 

0.28 

5 

-2 

9-12 

0.23 

15 

-2 

18-21 

0.23 

23-26 

0.28 

25 

-2 

28 

ground 

Note:  Depths  are  measured  down  from  the  pack  surface 


DATE. 

TIME. 


2/2$/ee 


ll:4S 


AMB  TEMP  (°C) 
SITE  NAME 


DEPTH  (ca) 

OENSin  (B/e«^ 

TE)«RATVJfle  t“C) 

1-4 

0.17 

5-8 

0.19 

5 

-7 

10-13 

0.17 

10 

-5 

14-17 

0.19 

15 

-3 

19 

yround 

DEPTH  (ca) 

CRYSTAL  SHE  (n) 

ICE  THICKNESS  (ca) 

0-17 

1  (tight) 

Note:  Depths  are  measured  down  from  the  pack  surface 
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UNDISTURBED  SNOW  CHARACTERIZATION 


riATF  2/2S/HB 

TIME _ liLiS _ 

AMB  TEMP  (°C)  _ 
SITE  NAME 


Note:  Depths  ere  meesureP  Poxn  from  the  peck  eurfece 

DATE  3/S/B8 _ _ 

TIME_iil£2 _ _ 

AMB  TEMP  (OC)  _ 

SITE  NAME  _ UC _ 


doth  (ca)  I  DEKSm  tj/ca^  I  TOreUTURE  (®C) 


0.21 


depth  lea)  I  CRYSTAL  SIZE  (aa)  I  ICE  THICKNESS  (caj 


1  (wet) 


1-1.5  (wet) 


Note:  Depths  are  measured  down  from  the  pscle  surface 
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UNDISTURBED  SNOW  CHARACTERIZATION 


HATF  3/3/aB _ 

TIME _ ii” _ 

AMB  TEMP  (Oc)  ^ 
SITE  NAME  EE _ 


DEP'm  Icb) 

CRYSTAL  SIZE  Ina) 

ICE  -miCKNESS  (ca) 

0 

4 

4-11 

1-1.5  (loose/drv) 

OEFm  ica) 

DENSITY  Ig/caT) 

TEMPERATURE  (“O 

0-3 

0.2G 

3 

-14 

4-7 

0.24 

7-10 

0.23 

9 

-S 

11 

ground 

Note:  Depths  are  measured  down  from  the  pack  surface 


hatf  3/3/ee  .. 
TIME 

AMB  TEMP  {°C) 
SITE  NAME 


DEPTH  Ice) 

CRYSTAL  SIZE  M 

ICE  TMICXNESS  (ca) 

0-3.5 

1  (tight/daap) 

3.5 

1  (snow) 

4.5-12 

1-1.5  (looaa/dBBi 

) 

DEPTH  (ca) 

DENSITT  (s/ea^ 

TEW>ERATURE  CO 

0-3 

0.15 

3-6 

0.27 

5 

0 

7-10 

0.20 

10 

-2 

12 

groufuS 

Note;  Depths  are  measured  down  from  the  pack  surface 
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UNDISTURBED  SNOW  CHARACTERIZATION 


nATF  3/4/B8 
TTMF  “:00 
AMB  TEMP  (°C) 
SITE  NAME 


DEPTH  (CB)  DENSITY  (g/cB^)  TEHPEHATUHE  (“0  DEPTH  lea)  CRYSTAL  SIZE  (aa)  ICE  THICKNESS  lea) 


DEPTH  (CB)  DENSITY  (g/ca^  TEM>EIUTURE  (®C)  DEPTH  (ca)  I  CRYSTAL  SIZE  ^  I  ICE  THICMCSS  (ca) 


Note:  Depths  ere  measured  down  from  the  peck  surface 
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UNDISTURBED  SNOW  CHARACTERIZATION 


63 


APPENDIX  E:  SHALLOW  SNOW  MOBILITY 
MODEL  CODE,  VERSION  2.0 
Includes  sample  output 

ie  I 

20  I  SSM/2.01 

30  I  (trig  tr  sinkage  and  snow  density  function  of  cont .  pr 

40  I  (pianual  input  of  all  vehicle  and  snou  data) 

50  I  (Liston’s  resis;  best-fit  power  curve  for  traction) 

60  I 

70  net=0 
80  RR=0 
90  DB=0 
100  go=0 
1 10  pl”0 
120  z=0 
130  sumz=0 

140  DISP  "vehicle’" 

150  INPUT  Ueh$ 

160  DISP  "wheeled  (u)  or  tracked  (t)’" 

170  INPUT  types 

180  DISP  "highest  station  ground  pressure  (kPa)*" 

190  INPUT  pnax 

200  pmax=pnax*0 . 1  !  convert  from  kPa  to  N/cm''2 
210  I 

220  !  ♦♦  establish  final  density  based  on  largest  footprint  pressure  •• 

230  sigmaf=0.5 

240  IF  pmax>21  THEN  sigmaf-0.55 
250  IF  pmax>35  THEN  sigmaf’'0.6 
260  IF  pmax>70  THEN  sigmaf*0.6S 
270  1 

280  DISP  "snou  depth  (cm)  »" 

290  INPUT  hO 

300  DISP  "snou  density  (g/cm''3)  •* 

310  INPUT  sigma0 
320  signal “Sigma0 
330  Sigma2*8igma0 
340  hl-h0 
350  I 

360  PRINT  'Uehicle:  "lUehS 
370  PRINT  ■  " 

380  PRINT  USING  formatl  i  "initial  state:  snow  depth  “'ih0i"  cm" 

390  PRINT  USING  format2  i  "  snou  density  “"isigmaGi"  g/cm"3 

400  I 

410  zmax»h0»( 1 -sigma0/sigmaf  )  !  calculate  maximum  sinkage 
420  IF  zmax<0  THEN  zmax“0 
430  I 

440  RAD  '  compute  in  radians  units  for  trig  functions 
450  I 
460  I 

470  I  ••  enter  vehicle  data  one  tire  station  at  a  time  •• 

480  DISP  “hou  many  wheel  or  track  stations  on  each  side  of  the  vehicle’" 
490  INPUT  N 
500  I 

510  FOR  I-l  TO  N 
520  DISP 

530  DISP  "station  ":I 

540  DISP  “wheel  radius  or  approx  radius  of  compacting  portion  of  track  cm) 
550  INPUT  rads 

560  DISP  “single  tire  or  track  width  at  this  location  (cm)  ■" 

570  INPUT  uid 

580  DISP  "single  tire  or  track  load  (N)  =" 

590  INPUT  loa 
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G00  OISP  "contact  area  <n‘2)  •“ 

610  INPUT  area 

620  IF  type$-"t"  THEN  GOTO  skip! 

630  DISP  "inflation  pressure  (kPa)  =" 

640  INPUT  pres 

650  pres*pres»0 . 1  '  convert  from  kPa  to  N/cn"2 
660  SOTO  skip2 

670  skipl:  pres* loa/area/ 1 0000  •  convert  from  N/fn“2  to  N/cn“2 
680  I 

690  skip2:  DISP  "driven’" 

700  INPUT  r$ 

710  IF  type$-"t"  THEN  GOTO  skip3 
720  OISP  "duals 
730  INPUT  R$ 

740  IF  R$<>"y"  AND  R$<>"Y"  THEN  SOTO  780 
750  uiid*uid»2 
780  loa*loa»2 
770  ! 

780  skip3:  DISP  "percent  of  width  compacting  virgin  snow  (?)  *" 

790  INPUT  prcnt 
800  I 

810  IF  pres>pl  AND  s igna2<s ignaf  AND  h0>0  THEN  GOTO  870  i  added  sink,  here’ 

820  z-0 

830  tempres'0 

840  IF  prcnt>0  AND  signa2>slgna0  THEN  GOTO  960 
850  resis=0 
860  GOTO  1010 

870  z*( pres/pmax  )"0.5»znax-sunz  !  calculate  additional  sinkage  this  station 

880  I 

890  I  set  rut  bottom  values 
900  h2-hl-z 

910  8igma2»s ignal /( 1 -z/ht  ) 

920  I. 

930  !  calculate  resistance  parameter  •• 

940  tempres"! /( sigma2-slgmal  )«L06< sigma2/siomal  )-l /sloma2  I  in  rut 
950  ! 

960  I  •*  calculate  motion  resistance  at  this  station  •• 

970  resi s'prcnt / 100»< pres»wld»h0«8iQma0  )  '  in  virgin  snow 

980  res i s*re8is»< 1 /< sigma2-8igma0  )»L06( signa2/sigma0 >- 1 /sigma2  )  Ivirgln  snow 
990  moreres” t empres»( 1 -prcnt / 1 00  )•< pres^wi d»h1 *8 Igma 1  )  '  in  rut 
1000  resis=resls+moreres 
1010  I 

1020  I  •*  calculate  gross  traction  at  this  station 
1030  trac*0.85l •( loa/area/1000 )'0.823  '  in  kPa 
1040  trac*trac»area» 1 000  !  in  N 
1050  I 

1060  IF  r$”"n"  OR  r$»"N“  THEN  trac“0  •  no  traction  if  not  driven 
1070  I 

1080  I  double  to  account  for  both  sides  of  the  vehicle  •• 

1090  arrea*2»area 
1 100  trac“2»trac 
1110  res  1 s”2»res i s 
1120  I 


1130 
1  140 
1150 
M60 
1170 
1180 
1  190 
1200 
1210 


••  print  output  for  this  station  •» 


PRINT  "  " 

PRINT  "station  "il 
PRINT  USING  format  I  i  " 
PRINT  USING  format  1  i  " 
PRINT  USING  formatl  i  “ 
PRINT  USING  formatl  i  " 

I 

I  ••  running  summation 


additional  sinkage  •".■z!” 
total  area*" tarrea I "  m“2" 
snow  resis  *"tresisi"  N" 
snow  trac  *"( tract"  N" 

for  whole  vehicle  •• 


cm 
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1220  RR*RR+resi5  I  sun  for  whole  vehicle 
1230  DB“DB+trac  '  sun  for  whole  vehicle 

1240  net-trac-resis  i  calculate  net  traction  for  individual  station 
1250  PRINT  USING  fornat 1  i  "  net  snow  traction  for  station  =  “ineti"  N“ 

1260  PRINT  USING  fornat 1  i  "  rut  bottom:  depth*"ih2i"  cm" 

1270  PRINT  USING  fornat2  i  "  dens i ty““ i s igna2  ;  ”  g/cn"3“ 

1280  go'go+net  i  sun  net  traction  for  vehicle 
1290  I 

1300  •  save  last  station  values  for  next  iteration  ♦♦ 

1310  IF  pres>pl  THEN  pl=pres 

1320  signal •5igna2 

1330  h1-h2 

1340  sunz=5un2+z 

1350  ! 

1360  NEXT  I 
1370  ! 

1380  I  ••  calculate  nobility  in  English  units  for  output 

1390  eRR»RR/4. 448222 

1400  eOB-DB/4. 448222 

1410  ego«go/4. 448222 

1420  ez"sunz/2.54 

1430  ! 

1440  I  *•  print  out  whole  vehicle  results  •• 

1450  PRINT  ■  " 

1460  PRINT  ■  ■ 

1470  PRINT  USING  fornatl  t  "  total  sinkage  for  vehicle”  "isunzi”  cn  (“tezi“  in 
)■ 

1480  PRINT  USING  fornatl  i  '  total  snown  resistance  •  “iRRi"  N  ("leRRi"  lb)“ 
1490  PRINT  USING  fornatl  i  ’  total  snow  traction  “’iDBt“  N  ("leDBi’  lb)” 

1500  PRINT  "  “ 

1510  PRINT  USING  fornatl  i  ”  net  traction  for  vehicle  ■  'iQoi'  N  ("legoi"  lb>" 
1520  PRINT  •  ■ 

1530  PRINT  •  “ 

1540  PRINT  •  '■ 

1550  I 

1560  fornatl:  IMAGE  3< K .DDOODD . DO  ) 

1570  fornat2:  IMAGE  3(K,0Z.0DD0) 

1580  I 
1590  END 
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Vehicle:  HMMUV 


initial  state:  snow  depth  •  5.00  cm 

snou  density  »  0.5500  g/cn"3 

station  I 

additional  sinkage  •  0.00  cm 

total  area=  0.1480  m“2 
snou  resis  •  0.00  N 

snou  trac  «  5487. G5  N 

net  snou  traction  for  station  «  5487.65  N 

rut  bottom:  depth*  0.00  cm 

density*  0.5500  3/cm''3 

station  2 

additional  sinkage  •  0.00  cm 

total  area*  0.1480  m"2 
snou  resis  *  0.00  N 

snou  trac  »  8844.84  N 

net  snow  traction  for  station  »  8844.84  N 

rut  bottom:  depth*  0.00  cm 

density*  0.5500  g/'cm"3 


total  sinkage  for  vehicle*  0.00  cm  (  0.00  in) 

total  snoun  resistance  *  0.00  N  (  0.00  lb) 

total  snow  traction  *  12312.49  N  (  2767.96  lb) 

net  traction  for  vehicle  *  12312.49  N  (  2767.96  lb) 
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Vehicle:  HEMTT  (24I/27G) 


initial  state:  snou  depth  »  25.00  cn 

snow  density  «  0.2200  q/ci*i''3 


station  1 

additional  sinkage  •=  14.03  cn 

total  area-  0.2980  n'2 
snou  resis  -  11515.55  N 
snou  tree  -  20514.00  N 
net  snou  traction  for  station  -  8998.45  N 

rut  bottom:  depth-  10.97  cm 

density-  0.5014  g/cn“3 

station  2 

additional  sinkage  -  0.00  cn 

total  area-  0.2980  ri"2 
snou  resis  -  0.00  N 

snou  trac  =  2078 1.93  N 
net  snou  traction  for  station  -  20781.93  N 

rut  bottom:  depth-  10.97  cm 

density-  0.5014  g/cm'3 


station  3 

additional  sinkage  ■  .97  cn 

total  area-  0.2980  m"2 
snow  resis  -  1230.94  N 

snow  trac  -  23084.34  N 
net  snou  traction  for  station  »  21853.40  N 

rut  bottom:  depth-  10.00  cn 

density-  0.5500  g/en'S 

station  4 

additional  sinkage  -  0.00  cm 

total  area-  0.2980  n''2 
snow  resis  -  0.00  N 

snou  trac  -  23214.12  N 
net  snou  traction  for  station  -  23214.12  N 

rut  botton:  depth-  10.00  cm 

density-  0.5500  g/cn*3 


total  sinkage  for  vehicle-  15.00  cm  (  5.91  In) 

total  snoun  resistance  -  12748.49  N  (  2885.52  lb) 

total  snou  traction  »  87574.39  N  (  19887.50  lb) 

net  traction  for  vehicle  -  74827.90  N  (  16821.98  lb) 
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Uehicle:  SUSU 

Initial  state:  snow  depth  “  9.00  cn 

snou  density  •  0.1480  g/cn“3 

station  t 

additional  sinkage  ”  6.30  cn 

total  area*  2.36  n'2 
snou  resis  ”  309.13  N 

snou  trac  “  16626.42  N 
net  snou  traction  for  station  »  16317.29  N 

rut  bottom:  depth*  2.70  cn 

density*  0.4930  q/cm*3 


station  2 

additional  sinkage  *  0.00  cm 

total  area*  2.36  m*2 
snou  resiS  •  0.00  N 

enow  trac  •  16527,49  N 
net  snou  traction  for  station  ■  16527.49  N 

rut  bottom:  depth*  2.70  cm 

density*  0.4930  o/cm*3 


total  sinkage  for  vehicle"  6.30  cm  (  2.48  in) 

total  snoun  resistance  *  309.13  N  <  69.49  lb) 

total  snow  traction  *  33153.90  N  (  7453.29  lb) 

net  traction  for  vehicle  *  32844.78  N  (  7383.80  lb) 


70 


REPORT  DOCUMENTATION  PAGE 


Form  Approved 
OMBNo.  0704-0188 


Public  repwriing  burden  for  this  collection  of  information  is  estimated  to  average  1  hour  per  response,  including  the  time  for  reviewing  instructions,  searching  existing  data  sources,  galhenng  and 
maintaining  the  data  needed,  and  completing  and  reviewing  the  collection  of  information  SerKJ  comments  regarding  this  burden  estimate  or  any  other  aspect  of  this  collection  of  information, 
Induding  suggestion  for  reducing  this  burden,  to  Washington  Headquarters  Services.  Directorate  for  Information  Operations  and  Reports.  12t5  Jefferson  Davjs  Highway,  Suite  1204.  Arlington. 
VA  22202*4302.  and  to  the  Office  of  Management  and  Budget.  Paperwork  Reduction  Project  (0704-0188).  Washington.  DC  20503. 


1.  AGENCY  USE  ONLY  (Leave  blank) 


4.  TITLE  AND  SUBTITLE 


2.  REPORT  DATE 

November  1990 


3.  REPORT  TYPE  AND  DATES  COVERED 


5.  FUNDING  NUMBERS 


Wheels  and  Tracks  in  Snow:  Validation  Study  of  the  CRREL  Shallow  Snow 
Mobility  Model 


6.  AUTHORS 

George  L.  Blaisdell.  Paul  W.  Richmond,  Sally  A.  Shoop,  Charles  E.  Green 
and  Russell  G.  Alser 


PE:  6.27.84A 
PR:  4A762784AT42 
TA:  CS 
WU:  040 


7.  PERFORMING  ORGANIZATION  NAME(S)  AND  ADDRESS(ES) 

U.S.  Army  Cold  Regions  Research  and  Engineering  Laboratory 
72  Lyme  Road 

Hanover,  New  Hampshire  03755-1290 


8.  PERFORMING  ORGANIZATION 
REPORT  NUMBER 


CRREL  Report  90-9 


9.  SPONSORING/MONITORING  AGENCY  NAME(S)  AND  ADDRESS(ES) 

Office  of  the  Chief  of  Engineers 
Washington,  DC  203 1 4- 1 000 


10  SPONSORING/MONITORING 
AGENCY  REPORT  NUMBER 


12a.  DISTRIBUTION// VAILABIUTY  STATEMENT 


Approved  for  public  release;  distribution  is  unlimited. 
Available  from  NTIS.  Springfield.  Virginia  22161 


13,  ABSTRACT  {Maximum  200  words) 

In  1986,  a  mobility  model  was  developed  for  predicting  the  traction  and  motion  resistance  of  both  wheeled  and 
tracked  vehicles  on  shallow  snow,  and  a  winter  field  season  was  dedicated  to  gathering  mobility  data  fora  diverse 
family  of  vehicles  (including  four  on  wheels  and  three  tracked)  to  validate  the  model.  The  original  version  of  the 
model,  SSMl  .0,  used  the  Mohr-Coulomb  shear  failure  equation  from  soil  mechanics  to  predict  gross  traction.  This 
required  input  of  the  snow  strength  parameters  c  and  <(>.  Motion  resistance  is  predicted  by  calculating  the  amount 
of  work  done  by  the  tire  in  compacting  snow  and  only  requires  snow  depth  and  density  values  as  input  snow 
properties.  Some  effort  was  expended  in  determining  an  easy  and  reliable  method  of  obtaining  snow  strength 
parameters.  The  model  was  originally  designed  to  use  an  initial  snow  density-snow  strength  relationship 
established  from  past  instrumented  vehicle  test  results.  Historically,  shear  annulus  apparati  have  been  used  to 
obtain  Mohr-Coulomb  strength  parameters.  A  comparison  of  snow  strength  obtained  via  these  three  methods 
(shear  annulus,  instrumented  vehicle,  calculated  from  initial  density  using  the  relationship  in  SSMl.O)  for 
individual  snow  covers  showed  no  agreement.  SSMl.O  assumed  that  snow  strength  parameters  for  mobility 
prediction  were  a  function  of  initial  snow  density;  however,  traction  is  developed  in  the  compacted  snow  under 
the  driving  element,  whose  strength  properties  bore  little  relation  to  those  of  the  initial  snow.  It  appears  that  the 
shear  strength  of  the  compacted  snow  is  essentially  a  constant  for  all  of  the  vehicles  and  snow  covers  tested  here. 


14.  SUBJECT  TERMS 

Computer  modeling 
Drawbar  pull 
Instrumented  vehicles 


17.  SECURITY  CLASSIFICATION 
OF  REPORT 

UNCLASSIFIED 


Mobility  predictions 
Motion  resistance 
Snow 


18.  SECURITY  CLASSIFICATION 
OF  THIS  PAGE 

UNCLASSIFIED 


Snow  strength 
Vehicle  traction 
Vehicles 


19.  SECURITY  CLASSIFICATION 
OF  ABSTRACT 

UNCLASSIFIED 


15.  NUMBER  OF  PAGES 

76 


16.  PRICE  CODE 


20.  LIMITATION  OF  ABSTRACT 


NSN  7540-01-280-5500 


Standard  Form  298  (Rev.  2-89) 

Pwscribwi  by  ANSI  Sid  Z39-18 
298-102 


13.  ABSTRACT  (confd). 

Based  on  this  finding,  a  new  traction  algorithm  was  developed,  resulting  in  the  creation  of  a  second  generation 
model,  SSM2.0.  This  algorithm  predicts  gross  traction,  on  the  average  for  the  vehicles  tested,  within  7%  of  the 
measured  value.  Motion  resistance  prediction  remains  unchanged  in  SSM2.0.  This  quantity  is  still  not  predicted 
with  a  desirable  level  of  accuracy. 


n 


'Pc  j .  .  r.f.vrf  JjMf  N'  f'BrjTpif.  jirj'i:  Hr:--'  .  2 


72 


